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Abstract  Every  three  years  the  IAU  Working  Group  on  Cartographic  Coordinates  and 
Rotational  Elements  revises  tables  giving  the  directions  of  the  poles  of  rotation  and  the 
prime  meridians  of  the  planets,  satellites,  minor  planets,  and  comets.  This  report  takes  into 
account  the  IAU  Working  Group  for  Planetary  System  Nomenclature  (WGPSN)  and  the  IAU 
Committee  on  Small  Body  Nomenclature  (CSBN)  definition  of  dwarf  planets,  introduces 
improved  values  for  the  pole  and  rotation  rate  of  Mercury,  returns  the  rotation  rate  of  Jupiter 
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to  a  previous  value,  introduces  improved  values  for  the  rotation  of  five  satellites  of  Saturn,  and 
adds  the  equatorial  radius  of  the  Sun  for  comparison.  It  also  adds  or  updates  size  and  shape 
information  for  the  Earth,  Mars’  satellites  Deimos  and  Phobos,  the  four  Galilean  satellites 
of  Jupiter,  and  22  satellites  of  Saturn.  Pole,  rotation,  and  size  information  has  been  added  for 
the  asteroids  (21)  Lutetia,  (511)  Davida,  and  (2867)  Steins.  Pole  and  rotation  information  has 
been  added  for  (2)  Pallas  and  (21)  Lutetia.  Pole  and  rotation  and  mean  radius  information 
has  been  added  for  (1)  Ceres.  Pole  information  has  been  updated  for  (4)  Vesta.  The  high 
precision  realization  for  the  pole  and  rotation  rate  of  the  Moon  is  updated.  Alternative  ori¬ 
entation  models  for  Mars,  Jupiter,  and  Saturn  are  noted.  The  Working  Group  also  reaffirms 
that  once  an  observable  feature  at  a  defined  longitude  is  chosen,  a  longitude  definition  origin 
should  not  change  except  under  unusual  circumstances.  It  is  also  noted  that  alternative  coor¬ 
dinate  systems  may  exist  for  various  (e.g.  dynamical)  purposes,  but  specific  cartographic 
coordinate  system  information  continues  to  be  recommended  for  each  body.  The  Working 
Group  elaborates  on  its  purpose,  and  also  announces  its  plans  to  occasionally  provide  limited 
updates  to  its  recommendations  via  its  website,  in  order  to  address  community  needs  for 
some  updates  more  often  than  every  3  years.  Brief  recommendations  are  also  made  to  the 
general  planetary  community  regarding  the  need  for  controlled  products,  and  improved  or 
consensus  rotation  models  for  Mars,  Jupiter,  and  Saturn. 

Keywords  Cartographic  coordinates  •  Longitude  •  Latitude  •  Rotation  axes  • 

Rotation  periods  •  Sizes  •  Shapes  •  Planets  •  Satellites  •  Dwarf  planets  •  Minor  planets  • 
Comets 


1  Introduction 

The  IAU  Working  Group  on  Cartographic  Coordinates  and  Rotational  Elements  of  the  Plan¬ 
ets  and  Satellites  was  established  as  a  consequence  of  resolutions  adopted  by  Commis¬ 
sions  4  and  16  at  the  IAU  General  Assembly  at  Grenoble  in  1976.  The  Working  Group 
became  a  joint  working  group  of  the  IAU  and  the  International  Association  of  Geodesy 
(IAG)  in  1985.  Due  to  a  lack  of  formal  communication  with  the  IAG  in  recent  years 
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that  affiliation  has  been  dropped  for  this  report.  It  may  be  re-established  in  the  future. 
Currently,  within  the  IAU,  the  Working  Group  is  a  joint  working  group  of  Divisions  I 
and  in,  and  not  part  of  any  commissions.  The  first  report  of  the  Working  Group  was 
presented  to  the  General  Assembly  at  Montreal  in  1979  and  published  in  the  Trans. 
IAU  17B,  72-79,  1980.  The  report  with  appendices  was  published  in  Celestial  Mechan¬ 
ics  22,  205-230,  1980.  The  guiding  principles  and  conventions  that  were  adopted  by  the 
Group  and  the  rationale  for  their  acceptance  were  presented  in  that  report  and  its  appen¬ 
dices.  The  second  report  of  the  Working  Group  was  published  in  the  Trans.  IAU  18B, 
151-162,  1983,  and  also  in  Celestial  Mechanics  29,  309-321,  1983.  In  2003  the  name 
of  the  Working  Group  was  shortened  to  the  Working  Group  on  Cartographic  Coordi¬ 
nates  and  Rotational  Elements.  The  following  table  summarizes  the  references  to  all  the 
reports. 


Report 

General  Assembly 

Celestial  Mechanics  and  Dynamical  Astronomy 

1 

Montreal  in  1979 

22,  205-230  (Davies  et  al.  1980) 

2 

Patras  in  1982 

29,  309-321  (Davies  et  al.  1983) 

3 

New  Delhi  in  1985 

39,  103-113  (Davies  et  al.  1986) 

4 

Baltimore  in  1988 

46,  187-204  (Davies  et  al.  1989) 

5 

Buenos  Aires  in  1991 

53,  377-397  (Davies  et  al.  1992) 

6 

The  Hague  in  1994 

63,  127-148  (Davies  et  al.  1996) 

7 

Kyoto  in  1997 

No  report 

8 

Manchester  in  2000 

82,  83-110  (Seidelmann  et  al.  2002) 

9 

Sydney  in  2003 

91,  203-215  (Seidelmann  et  al.  2005) 

10 

Prague  in  2006 

98,  155-180  (Seidelmann  et  al.  2007) 

11 

Rio  de  Janeiro  in  2009 

This  paper 

Reprints  and  preprints  of  the  previous  reports  and  this  report  can  be  found  at  the  Working 
Group  web  site:  http://astrogeology.usgs.gov/Projects/WGCCRE.  Previous  reports  are  also 
available  at  the  web  site:  http://www.springerlink.com/content/100246. 

The  original  impetus  for  the  Working  Group  was  stated  in  an  IAU  Resolution:  “to  avoid  a 
proliferation  of  inconsistent  cartographic  and  rotational  systems,  there  is  a  need  to  define  the 
cartographic  and  rotational  elements  of  the  planets  and  satellites  on  a  systematic  basis  and  to 
relate  the  new  cartographic  coordinates  rigorously  to  the  rotational  elements”  (International 
Astronomical  Union  (IAU)  1977,  p.  144).  Since  its  first  report  (Davies  et  al.  1980),  this 
Working  Group  has  addressed  this  need  and  its  purpose  has  remained  essentially  unchanged, 
except  for  the  recognition  of  the  need  to  address  the  same  issues  for  small  bodies  of  the  Solar 
System,  beginning  with  the  2003  report. 

Therefore  in  actual  execution,  the  Working  Group  sees  its  mission  as  this:  making  rec¬ 
ommendations  that  define  and  relate  the  coordinate  systems  of  Solar  System  bodies  to  their 
rotational  elements  to  support  making  cartographic  products  (i.e.  “mapping”)  of  such  bod¬ 
ies.  The  working  group  incorporates  any  reasonable  and  peer  reviewed  improved  determina¬ 
tions  that  follow  previously  established  conventions,  or  possibly  in  some  cases  may  decide 
between  different  such  determinations.  Because  of  the  lack  of  the  necessary  resources,  how¬ 
ever,  the  Working  Group  does  not  verify  or  validate  such  determinations.  Our  recommen¬ 
dations  are  from  the  Working  Group  alone  and,  if  only  for  reasons  of  practicality,  are  not 
recommendations  from  the  full  IAU.  The  Working  Group  has  no  “enforcement”  mecha¬ 
nism  to  assure  that  its  recommendations  are  followed — the  value  of  these  recommendations 
is  only  from  their  development  by  international  consensus  and  adoption  by  the  planetary 
community. 
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The  main  use  of  these  recommendations  should  be  in  cases  where  standardization  is  use¬ 
ful.  It  is,  of  course,  not  our  intention  to  limit  science  or  the  state  of  the  art.  If  for  cartographic 
or  other  purposes  these  recommendations  can  be  followed  without  affecting  the  quality  of 
the  product,  then  they  should  be  followed  so  that  products  can  be  more  easily  compared 
and  multiple  datasets  appropriately  registered.  If  a  user  has  sufficient  data  to  update  the  rec¬ 
ommended  models  or  values  used  so  that  any  cartographic  product  would  be  significantly 
improved,  then  obviously — following  the  conventions  described  here — such  updates  should 
be  used.  (It  may  also  be  useful  to  make  alternative  products  using  old  and  new  models  for 
comparison  purposes). 

Because  there  will  inevitably  be  some  delay  before  our  next  triennial  report  is  published 
(or  our  website  updated — see  below),  this  type  of  action  is  almost  always  necessary  when 
updated  parameter  estimates  are  derived.  We  encourage  the  publication  of  such  updates  in  the 
peer-reviewed  literature  at  the  earliest  opportunity.  This  is  desirable  both  so  that  the  Working 
Group  can  consider  them  and  as  appropriate  recommend  their  use  and  properly  reference 
them  in  our  next  report,  but  also  so  others  will  become  aware  of  such  updates,  and  can  use 
them  if  necessary  in  the  interim. 

At  the  request  of  various  individuals  and  missions,  in  the  coming  triennium  the  Working 
Group  will  consider  providing  limited  updates  to  its  recommendations  via  its  website.  This 
is  in  order  to  address  a  cited  need  to  update  recommendations  more  often  than  every  3  years, 
e.g.  for  use  by  operational  missions  and  for  updated  cartographic  products.  Details  of  this 
procedure  are  still  to  be  worked  out.  However,  our  tentative  plan  is  to  determine  approxi¬ 
mately  every  6  months  whether  such  updates  are  necessary  and,  if  so,  announce  them  on 
our  web  site.  We  also  plan  to  begin  listing  on  the  website — without  any  recommendations 
and  for  information  purposes  only — new  recently  published  and  (preferably)  peer-reviewed 
determinations  related  to  Solar  System  coordinate  systems.  Note  that  the  posting  of  recom¬ 
mendations  to  our  web  site  is  not  intended  to  supersede  the  need  for  our  triennial  reports. 
We  plan  to  reserve  the  bulk  of  any  new  recommendations  or  changes  to  our  recommenda¬ 
tions  for  these  reports,  only  placing  time-critical  ones  on  our  website.  In  our  next  report,  we 
will  consider  the  usefulness  of  these  procedures  and  whether  they  should  continue.  Input  for 
such  updates  (whether  for  Working  Group  consideration  or  information  only)  and  comment 
on  these  procedures  from  the  community  is  welcome,  particularly  from  missions  and  other 
working  groups. 

The  2003  report  introduced  and  recommended  a  consistent  system  of  coordinates  for  both 
minor  planets  (asteroids)  and  comets,  and  in  this  report  we  extend  it  to  cover  dwarf  planets. 
This  system  is  not  the  same  as  the  system  for  planets  and  satellites,  which  was  not  changed. 
It  is  recognized  that  the  existence  of  two  different  systems  has  the  potential  for  confusion, 
but  the  methods  required  for  dwarf  planets,  minor  planets,  their  satellites,  and  comets  differ 
sufficiently  to  justify  the  use  of  two  different  systems.  This  report  includes  descriptions  of 
the  two  systems;  one  for  planets  and  satellites  (Sects.  2,3,4,  6)  and  another  for  dwarf  planets, 
minor  planets,  their  satellites,  and  comets  (Sects.  5  and  7).  Rotational  elements  (body  orien¬ 
tation  in  inertial  space)  are  covered  first  (Sects.  2-5),  and  then  cartographic  coordinates,  e.g. 
latitude,  longitude,  and  body  shape  (Sects.  6,  7).  Brief  recommendations  from  the  Working 
Group  complete  this  report  (Sect.  8).  For  the  purpose  of  assigning  bodies  such  as  Pluto  and 
Ceres  to  these  tables,  this  report  assumes  that  dwarf  planets  are  the  primary  bodies  on  the 
list  maintained  by  the  IAU  Working  Group  for  Planetary  System  Nomenclature  (WGPSN) 
and  the  IAU  Committee  on  Small  Body  Nomenclature  (CSBN)  (2010).  The  Appendix  lists 
changes  made  since  the  previous  report. 
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2  Definition  of  Rotational  Elements  for  Planets  and  Satellites 

Planetary  coordinate  systems  are  defined  relative  to  their  mean  axis  of  rotation  and  var¬ 
ious  definitions  of  longitude  depending  on  the  body.  The  longitude  systems  of  most  of 
those  bodies  with  observable  rigid  surfaces  have  been  defined  by  references  to  a  sur¬ 
face  feature  such  as  a  crater.  Approximate  expressions  for  these  rotational  elements  with 
respect  to  the  International  Celestial  Reference  Frame  (ICRF)  (Ma  et  al.  1998)  have  been 
derived.  The  ICRF  is  the  reference  frame  of  the  International  Celestial  Reference  Sys¬ 
tem  and  is  itself  epochless.  There  is  a  small  (well  under  0. 1  arcsecond)  rotation  between 
the  ICRF  and  the  mean  dynamical  frame  of  J2000.0.  The  epoch  J2000.0,  which  is  JD 
2451545.0  (2000  January  1  12  hours),  TDB,  is  the  epoch  for  variable  quantities,  which  are 
expressed  in  units  of  days  (86400  SI  seconds)  or  Julian  centuries  of  36525  days.  Barycen- 
tric  Dynamical  Time  (TDB)  is  the  reference  time  scale  for  time  dependent  variables.  TDB 
was  clarified  in  definition  at  the  IAU  General  Assembly  of  2006  in  Prague.  TDB,  some¬ 
times  called  Teph,  is  roughly  equivalent  to  Terrestrial  Time  (TT)  in  epoch  and  rate.  UTC, 
TCB,  and  TCG  differ  from  TT  in  epoch  and  rate.  For  more  information  on  reference  sys¬ 
tems  and  time  scales  see  Kovalevsky  and  Seidelmann  (2004),  http://www.iers.org,  http:// 
rorf.usno.navy.mil/ICRF/,  or  http://www.usno.navy.mil/USNO/astronomical-applications/ 
astronomical-information-center/icrs-narrative. 

The  north  pole  is  that  pole  of  rotation  that  lies  on  the  north  side  of  the  invariable  plane  of 
the  Solar  System.  The  direction  of  the  north  pole  is  specified  by  the  value  of  its  right  ascension 
ao  and  declination  So.  With  the  pole  so  specified,  the  two  intersection  points  of  the  body’s 
equator  and  the  ICRF  equator  are  ao  ±  90°.  We  choose  one  of  these,  ao  +  90°,  and  define 
it  as  the  node  Q.  Suppose  the  prime  meridian  has  been  chosen  so  that  it  crosses  the  body’s 
equator  at  the  point  B.  We  then  specify  the  location  of  the  prime  meridian  by  providing  a 
value  for  W,  the  angle  measured  easterly  along  the  body’s  equator  between  the  node  Q  and 
the  point  B  (see  Fig.  1 ).  The  right  ascension  of  the  point  Q  is  90°  +  ao  and  the  inclination  of 
the  planet’s  equator  to  the  celestial  equator  is  90°  —  So-  As  long  as  the  planet,  and  hence  its 
prime  meridian,  rotates  uniformly,  W  varies  nearly  linearly  with  time.  In  addition,  ao.So,  and 
W  may  vary  with  time  due  to  a  precession  of  the  axis  of  rotation  of  the  planet  (or  satellite). 


Fig.  1  Reference  system  used  to 
define  orientation  of  the  planets 
and  their  satellites 
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If  W  increases  with  time,  the  planet  has  a  direct  (or  prograde)  rotation,  and,  if  W  decreases 
with  time,  the  rotation  is  said  to  be  retrograde. 

In  the  absence  of  other  information,  since  most  satellites  fall  into  this  category,  the  axis  of 
rotation  is  assumed  to  be  normal  to  the  mean  orbital  plane  of  the  planet  or  the  satellite.  For 
many  of  the  satellites,  it  is  assumed  that  the  rotation  rate  is  equal  to  the  mean  orbital  period 
(i.e.  synchronous  rotation),  but  in  some  cases  such  an  assumption  still  needs  to  be  validated. 

The  angle  W  specifies  the  ephemeris  position  of  the  prime  meridian.  For  planets  or  sat¬ 
ellites  without  any  accurately  observable  fixed  surface  features,  the  adopted  expression  for 
W  defines  the  prime  meridian  and  is  not  subject  to  correction  for  this  reason.  However,  the 
rotation  rate  may  be  redefined  for  other  reasons.  Where  possible,  however,  the  cartographic 
position  of  the  prime  meridian  is  defined  by  a  suitable  observable  feature,  and  so  the  con¬ 
stants  in  the  expression  W  =  Wo  +  Wd,  where  d  is  the  interval  in  days  from  the  standard 
epoch,  are  chosen  so  that  the  ephemeris  position  follows  the  motion  of  the  cartographic 
position  as  closely  as  possible;  in  these  cases  the  expression  for  W  may  require  emendation 
in  the  future.  When  new  higher  accuracy  mapping  is  done,  the  longitude  of  the  fixed  feature 
should  be  maintained  and  a  new  value  for  Wo  derived,  the  results  published  in  peer-reviewed 
literature,  and  the  result  reported  to  this  Working  Group  for  possible  adoption.  For  bodies 
where  they  are  in  use,  longitude  defining  features  are  noted  in  the  footnotes  to  Tables  1,  2, 
and  3. 

The  Working  Group  would  like  to  emphasize — as  it  did  in  the  introduction  to  its  first  report 
(Davies  et  al.  1980,  p.  73) — that  once  an  observable  feature  at  a  defined  longitude  is  chosen, 
the  longitude  definition  origin  should  not  change  except  under  unusual  circumstances  (such 
as  perhaps  a  change  in  or  loss  of  the  feature).  This  implies  that  once  such  a  feature  has  been 
adopted,  a  return  to  a  value  of  Wo  defined  by  some  other  method  (e.g.  the  principal  axes  of 
inertia  for  resonantly  or  synchronously  rotating  bodies  such  as  Mercury  (Margot  2009),  or 
lovian  or  Saturnian  satellites)  should  be  avoided.  Note,  however,  that  this  does  not  preclude 
the  use  of  smaller  or  more  precisely  determined  features,  multiple  features,  or  even  human 
artifacts  to  define  longitude — as  long  as  the  original  definition  is  maintained  to  within  the 
accuracy  of  previous  determinations.  An  example  is  the  redefinition  of  the  origin  for  longi¬ 
tude  for  Mars  from  the  large  feature  then  known  as  Sinus  Meridiani  to  the  small  crater  Airy-0 
(de  Vaucouleurs  et  al.  1973). 

Recommended  values  of  the  constants  in  the  expressions  for  ao,  So,  and  W,  in  celestial 
equatorial  coordinates,  are  given  for  the  planets  and  satellites  in  Tables  1  and  2.  In  general, 
these  expressions  should  be  accurate  to  one-tenth  of  a  degree;  however,  two  decimal  places 
are  given  to  assure  consistency  when  changing  coordinates  systems.  Zeros  have  sometimes 
been  added  to  rate  values  (W)  for  computational  consistency  and  are  not  an  indication  of 
significant  accuracy.  Additional  decimal  places  are  given  in  the  expressions  for  Mercury,  the 
Moon,  Mars,  Saturn,  and  Uranus,  reflecting  the  greater  confidence  in  their  accuracy.  Expres¬ 
sions  for  the  Sun  and  Earth  are  given  to  a  similar  precision  as  those  of  the  other  bodies  of  the 
Solar  System  and  are  for  comparative  purposes  only. 

These  recommendations  are  not  intended  to  imply  that  other  coordinate  systems  with  dif¬ 
ferent  rotational  elements  should  not  be  used  for  planetary  bodies  for  other  than  cartographic 
or  provisional  purposes.  For  example  it  is  recognized  that  the  use  of  dynamical  coordinate 
systems  such  as  those  tied  to  a  body’s  principal  axis  may  be  needed  for  computational  pur¬ 
poses  or  for  important  dynamical  work.  The  use  of  such  coordinate  systems  is  common, 
for  example  as  described  immediately  below  for  the  Moon,  and  also  for  Mercury  (Margot 
2009).  It  is  also  possible,  depending  for  example  on  the  observational  mode  and  accuracy, 
that  a  body  fixed  coordinate  frame  can  at  times  be  defined  relative  to  inertial  space  at  a  higher 
level  of  accuracy  than  to  a  surface  feature  fixed  frame.  Mercury  is  again  an  example,  where 
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Table  1  Recommended  values  for  the  direction  of  the  north  pole  of  rotation  and  the  prime  meridian  of  the 
Sun  and  planets 

ao,  <5o  are  ICRF  equatorial  coordinates  at  epoch  J2000.0. 

Approximate  coordinates  of  the  north  pole  of  the  invariable  plane  are  o?o  = 
T  =  interval  in  Julian  centuries  (of  36525  days)  from  the  standard  epoch 
d  =  interval  in  days  from  the  standard  epoch 

The  standard  epoch  is  JD  2451545.0,  i.e.  2000  January  1  12 hours  TDB 

273°. 85.  i0  =  66°. 99 

Sun 

<*o  =  286°. 13 

S0  =  63°. 87 

W  =  84°. 176+  14°. 1844000 d 

(a) 

Mercury 

a0  =  281.0097-0.03287 

«50  =  61.4143  -0.00497- 

W  =  329.5469  +6.1385025(7 

+0°. 00993822  sin  (Ml) 

—0°. 00104581  sin  (M2) 

-0°. 00010280  sin  (M3) 

-0°. 00002364  sin (M4) 

-0°. 00000532  sin  (M5) 

where 

Ml  =  174°. 791086  +  4°  .092335(7 

M2  =  349°  .582171  +  8°. 184670(7 

M3  =  164°. 373257  +  12°.277005(7 

M4  =  339°.  164343  +  16°.369340(7 

M5  =  15  3°. 95 5429  +  20+461675(7 

(b) 

Venus 

a0  =  272.76 

S0  =  67.16 

W  =  160.20-  1.4813688(7 

(c) 

Earth 

o-o  =  0.00 -0.641 7- 
£>0  =  90.00  -  0.5577 

W  =  190.147  +  360.9856235£/ 

Mars 

ccq  =  317.68143  —  0.10617 

So  =  52.88650  -  0.06097 

W  =  176.630  +  350.89198226(7 

(d) 

Jupiter 

ff0  =  268.056595  -  0.0064997  +  0° .0001 17sin  Ja  +  0° .000938  sin  Jb 

+  0.001432  sin  Jc  +  0.000030  sin  Jd  +  0.002150  sin  Je 
«0  =  64.495303  +  0.0024137  +  0.000050  cos  Ja  +  0.000404  cos  Jb 

+  0.000617  cos  Jc  —  0.000013  cos  Jd  +  0.000926  cos  Je 

W  =  284.95  +  870.5360000(7 

(e) 

where  Ja  =  99°.360714  +  4850+40467.  Jb  =  175°.895369  +  1 191°.96057, 
Jc  =  300°.323162  + 262+54757,  Jd  =  1 14°.012305  +  6070°.24767, 
Je  =  49°.511251  +64°.30007 
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Table  1  continued 

Saturn 

cto  =  40.589  -  0.0367" 

8q  =  83.537  -  0.004 T 

IV  =  38.90  +  810.7939024c/ 

Uranus 

ag  =  257.311 

S0  =  -15.175 

(e) 

W  =  203.81  -  501.1600928J 

Neptune 

ofQ  =  299.36  +  0.70  sin  N 

<50  =  43.46  -0.51  cos  N 

(e) 

IV  =  253.18  +  536.3128492J  -  0.48  sin  N 

N  =  357.85  + 52.3 16T 

(e) 

(a)  The  equation  W  for  the  Sun  is  now  corrected  for  light  travel  time  and  removing  the  aberration  correction. 
See  the  Appendix  in  Seidelmann  et  al.  (2007) 

(b)  The  20°  meridian  is  defined  by  the  crater  Hun  Kal 

(c)  The  0°  meridian  is  defined  by  the  central  peak  in  the  crater  Ariadne 

(d)  The  0°  meridian  is  defined  by  the  crater  Airy-0 

(e)  The  equations  for  Wfor  Jupiter,  Saturn,  Uranus  and  Neptune  refer  to  the  rotation  of  their  magnetic  fields 
(System  III).  On  Jupiter,  System  I  ( Wj  =  67°.  1  +  877°.900J)  refers  to  the  mean  atmospheric  equatorial 
rotation;  System  II  (Wjj  =  43°. 3  +  870°.270<i)  refers  to  the  mean  atmospheric  rotation  north  of  the  south 
component  of  the  north  equatorial  belt,  and  south  of  the  north  component  of  the  south  equatorial  belt 


Table  2  Recommended  values  for  the  direction  of  the  north  pole  of  rotation  and  the  prime  meridian  of  the 
satellites 

ao,  <$o>  T,  and  d  have  the  same  meanings  as  in  Table  1  (epoch  JD  2451545.0,  i.e.  2000  January  1  12hours  TDB) 

Earth: 

Moon  a0  =  269°. 9949  +  0° .003 1 T  -  3°. 8787  sin  £1  -  0°.1204  sin  £2  (a) 

+  0.0700  sin  £3  -  0.0172  sin  £4  +  0.0072  sin  £6 
-0.0052  sin  £10+0.0043  sin  £13 

<$0  =  66.5392  +  0.0130T  +  1.5419  cos  £1  +  0.0239  cos  £2 

—  0.0278  cos  £3  +  0.0068  cos  £4  —  0.0029  cos  £6 
+  0.0009  cos  £7  +  0.0008  cos  £10  —  0.0009  cos  £13 

W  =  38.3213  +  13.17635815^  -  1.4x  I0~nd2  +  3.5610  sin  £1 
+  0.1208  sin  £2  -  0.0642  sin  £3  +  0.0158  sin  £4 
+  0.0252  sin  £5  -  0.0066  sin  £6  -  0.0047  sin  £7 

—  0.0046  sin  £8  +  0.0028  sin  £9  +  0.0052  sin  £  10 

+  0.0040  sin  £11 +0.0019  sin  £12  -  0.0044  sin  £13 
where 

£1  =  125°. 045  -  0°.0529921d,  £2  =  250° .089  -  0°.1059842</,  £3  =  260° .008  +  13°.0120009d, 

£4=  176.625  +  13.3407154J,  £5  =  357.529  + 0.9856003d,  £6=  31 1.589  + 26.4057084d, 

£7=  134.963+  13.0649930J,  £8  =  276.617  +  0.3287146*/,  £9=34.226+  1.7484877 d, 

£  10  =  15. 134  -  0. 1589763d,  £11  =  1 19.743  +  0.0036096d,  £  12  =  239.961  +  0. 1643573 d, 

£  13  =  25.053  +  12.9590088^/ 
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Table  2  continued 


Mars: 


I 

Phobos 

ct0  =  317.68  -  0.1087  +  1.79  sin  M 1 

S0  =  52.90-0.0617-  1.08  cos  Ml 

W  =  35.06  +  1 128.84458507  +  8.86472 
-  1.42  sin  Ml -0.78  sin  M2 

II 

Deimos 

ct0  =  316.65  -  0.1087  +  2.98  sin  M3 

S0  =  53.52  -  0.0617  -  1.78  cos  M3 

W  =  79.41  +  285.1618970(7  —  0.52072 
—  2.58  sin  M3  +  0.19  cos  M3 

where  Ml 

=  169°. 51  —  0°.43576407,  M2  =  192°  .93+  1128°. 4096700c?  +  8°.86472, 

M3 

=  53°. 47  —  0°. 0181510c? 

Jupiter: 

XVI 

Metis 

a0  =  268.05  -  0.0097 

S0  =  64.49  +  0.0037 

W  =  346.09  +  1221.2547301c? 

XV 

Adrastea 

ffO  =  268.05  -  0.0097 

S0  =  64.49  +  0.0037 

W  =  33.29  +  1206.9986602(7 

V 

Amalthea 

ot0  =  268.05  -0.0097  -  0.84  sin  71+0.01  sin  271 

S0  =  64.49  +  0.0037  -  0.36  cos  7 1 

W  =  231.67  +  722.63145607  +  0.76  sin  71  -  0.01  sin  271 

XIV 

Thebe 

ffO  =  268.05  -  0.0097  -  2.11  sin  72  +  0.04  sin  272 

S0  =  64.49  +  0.0037  -  0.91  cos  72  +  0.01  cos  272 

W  =  8.56  +  533.70041007  +  1.91  sin  72  -  0.04  sin  272 

I 

Io 

a0  =  268.05  -  0.0097  +  0.094  sin  73  +  0.024  sin  74 

S0  =  64.50  +  0.0037  +  0.040  cos  73  +  0.01 1  cos  74 

W  =  200.39  +  203.48895387  -  0.085  sin  73  -  0.022  sin  74 

II 

Europa 

a0  =  268.08  -  0.0097  +  1.086  sin  74  +  0.060  sin  75 
+  0.015  sin  76  +  0.009  sin  77 

S0  =  64.51  +0.0037 +  0.468  cos  74 +  0.026 cos  75 
+  0.007  cos  76  +  0.002  cos  77 

W  =  36.022  +101.37472357  -  0.980  sin  74  -  0.054  sin  75 
—  0.014  sin  76  —  0.008  sin  77 

III 

Ganymede 

a0  =  268.20  -  0.0097  -  0.037  sin  74  +  0.431  sin  75  +  0.091  sin  76 

S0  =  64.57  +  0.0037  -  0.016 cos  74  +  0.186 cos  75  +  0.039  cos  76 

W  =  44.064  +  50.31760817  +  0.033  sin  74  -  0.389  sin  75  -  0.082  sin  76 

IV 

Callisto 

a0  =  268.72  -  0.0097  -  0.068  sin  75  +  0.590  sin  76  +  0.010  sin  78 

So  =  64.83  +  0.003T  -  0.029  cos  75  +  0.254  cos  76  -  0.004  cos  78 

W  =  259.51  +  21.57107157  +  0.061  sin  75  -  0.533  sin  76  -  0.009  sin  78 

where  71  =  73°.32  +  91472° .97.  72  =  24°.62  +  45137°.27, 

73  =  283° .90  +  4850° .77,  74  =  355.80+  1191.37,  75  =  119.90  +  262.17, 
76  =  229.80  +  64.37.  77  =  352.25  +  2382.67,  78  =  1 13.35  +  6070.07 

Saturn: 

XVIII 

Pan 

<*0  =  40.6  -  0.0367 

S0  =  83.5  -  0.0047 

W  =  48.8  +  626.04400007 

XV 

Atlas 

ct0  =  40.58  -  0.0367 

S0  =  83.53  -  0.0047 

W  =  137.88  +  598.30600007 

(b) 


(c) 


(d) 


(e) 
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Table  2 

continued 

XVI 

Prometheus 

no  =  40.58  -  0.036 7 

S0  =  83.53  -  0.0047 

W  =  296.14  +  587.289000*7 

XVII 

Pandora 

n0  =  40.58  -  0.036 7 

S0  =  83.53  -  0.0047 

W  =  162.92  +  572.7891000rf 

XI 

Epimetheus 

ao  =  40.58  -  0.0367  -  3.153  sin  SI  +  0.086  sin  2S1 

S0  =  83.52  -  0.0047  -  0.356  cos  SI  +  0.005  cos  2S1 

W  =  293.87  +  518.4907239(7  +  3.133  sin  SI  -  0.086  sin  2S1 

(f) 

X 

Janus 

a0  =  40.58  -  0.0367  -  1.623  sin  S2  +  0.023  sin  2S2 

S0  =  83.52  -  0.0047  -  0.183  cos  S2  +  0.001  cos  2S2 

W  =  58.83  +  518.2359876(7  +  1.613  sin  S2  -  0.023  sin  2S2 

(f) 

I 

Mimas 

o-O  =  40.66  -  0.0367  +  13.56  sin  S3 

So  =  83.52  -  0.0047  -  1.53  cos  S3 

W  =  333.46  +  381.9945550*7  -  13.48  sin  S3  -  44.85  sin  S5 

(g) 

II 

Enceladus 

n0  =  40.66  -  0.0367 

S0  =  83.52  -  0.0047 

W  =  6.32  +  262.7318996*7 

(h) 

III 

Tethys 

no  =  40.66  —  0.0367  +  9.66  sin  S4 

So  =  83.52  -  0.0047  -  1 .09  cos  S4 

W  =  8.95  +  190.6979085(7  -  9.60  sin  S4  +  2.23  sin  S5 

(0 

XIII 

Telesto 

no  =  50.51  -  0.0367 

S0  =  84.06  -  0.0047 

W  =  56.88  +  190.6979332*7 

(f) 

XIV 

Calypso 

n0  =  36.41  -  0.0367 

S0  =  85.04  -  0.0047 

W  =  153.51  +  190.6742373*7 

(f) 

IV 

Dione 

n0  =  40.66  -  0.0367 

S0  =  83.52  -  0.0047 

W  =  357.6+  131.5349316*7 

(j) 

XII 

Helene 

n0  =  40.85  -  0.0367 

S0  =  83.34  -  0.0047 

W  =  245.12+  131.6174056*7 

V 

Rhea 

no  =  40.38  —  0.0367  +  3.10  sin  S6 

So  =  83.55  -  0.0047  -  0.35  cos  S6 

W  =  235.16  +  79.6900478*7  -  3.08  sin  S6 

00 

VI 

Titan 

n0  =  39.4827 

S0  =  83.4279 

W  =  186.5855  +  22.5769768*7 
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Table  2 

continued 

VIII 

Iapetus 

cro  =  318.16  —  3.9497’ 

S0  =  75.03  -  1.1437’ 

W  =  355.2  +  4.5379572<7  (1) 

IX 

Phoebe 

cto  =  356.90 
<50  =  77.80 

W  =  178.58  +  931.639(7 

where  S'!  = 

353°. 32  +  75706°.7 T,  52  =  28°.72  +  75706°.77,  53  =  177° .40  -  36505°.57 

Uranus: 

S4  =  300.00  -  7225.97/  55  =  316.45  +  506.27’,  S6  =  345.20  -  1016.37’ 

VI 

Cordelia 

a o  =  257.31  -0.15  sin  C/1 

S0  =  -15.18 +  0.14  cos  £71 

W  =  127.69  -  1074.5205730(7  -  0.04  sin  U 1 

VII 

Ophelia 

ao  =  257.31  -0.09  sin  C/2 

S0  =  -15.18  +  0.09  cos  C/2 

W  =  130.35  -  956.4068150(7  -  0.03  sin  £72 

VIII 

Bianca 

ao  =  257.31  -0.16  sin  C/3 
<50  =  -15.18 +  0.16  cos  C/3 

W  =  105.46  -  828.3914760(7  -  0.04  sin  C/3 

IX 

Cressida 

ck0  =  257.31  -0.04  sin  £74 

S0  =  -15.18 +  0.04  cos  C/4 

W  =  59.16  -  776.5816320(7  -  0.01  sin  C/4 

X 

Desdemona 

cto  =  257.31  -0.17  sin  C/5 

30  =  -15.18 +  0.16  cos  C/5 

W  =  95.08  -  760.0531690c/  -  0.04  sin  C/5 

XI 

Juliet 

a0  =  257.31  -0.06  sin  C/6 

<50  =  -15.18  +  0.06  cos  C/6 

W  =  302.56  -  730.1253660c/  -  0.02  sin  U6 

XII 

Portia 

a0  =  257.31  -  0.09  sin  C/7 

S0  =  -15.18  +  0.09  cos  C/7 

W  =  25.03  -  701.4865870c/  -  0.02  sin  C/7 

XIII 

Rosalind 

cc0  =  257.31  -0.29  sin  C/8 

S0  =  -15.18  +  0.28  cos  U8 

W  =  314.90  -  644.6311260c/  -  0.08  sin  US 

XIV 

Belinda 

cto  =  257.31  -0.03  sin  C/9 
<50  =  -15.18  +  0.03  cos  C/9 

W  =  297.46  -  577.3628170c/  -  0.01  sin  C/9 

XV 

Puck 

ck0  =  257.31  -0.33  sin  U 10 
<50  =  -15.18  +  0.31  cos  C/10 

W  =  91.24  -  472.5450690c/  -  0.09  sin  C/10 

V 

Miranda 

CK0  =  257.43  +  4.41  sin  C/ll  -  0.04  sin  2C/11 

S0  =  -15.08  +  4.25  cos  C/11-0.02  cos2£/ll 

W  =  30.70  -  254.6906892c/  -  1.27  sin  C/12  +  0.15  sin  21/12 

+  1.15  sin  U 11  -  0.09  sin  2C/1 1 
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Table  2  continued 


I 


II 


III 


IV 


Neptune 

III 


IV 


V 


VI 


VII 


VIII 


I 


Ariel  ag  =  257.43  +  0.29  sin  C/13 
<50  =  -15.10  +  0.28  cos  U 1 3 

IV  =  156.22  -  142.8356681(1  +  0.05  sin  1712  +  0.08  sin  1713 
Umbriel  ag  =  257.43  +  0.21  sin  1714 
<50  =  -15.10  +  0.20  cos  1714 

IV  =  108.05  -  86.8688923<7  -  0.09  sin  1712  +  0.06  sin  1714 
Titania  ag  =  257.43  +  0.29  sin  17 15 
<50  =  -15.10  +  0.28  cost/15 
IV  =  77.74  -  41.3514316(7  +  0.08  sin  1/15 
Oberon  ag  =  257.43  +  0.16  sin  1716 
Sg  =  -15.10 +  0.16  cos  1/16 
IV  =  6.77  -  26.7394932 d  +  0.04  sin  17 16 

where  1/1  =  1 15°.75  +  54991°.877\  172  =  141°,69  +  41887°.667,  1/3=  135o.03+29927°.357, 
174  =  61.77  +  25733.5977  t/5  =  249.32  +  24471.4677  176  =  43.86  +  22278.417. 

U1  =  77.66  +  20289.427  t/8  =  157.36  +  16652.767,  1/9  =  101.81  +  12872.637, 

1/10  =  138.64  +  8061.817,  1/11  =  102.23-  2024.227.  1712  =  316.41  +2863.967, 

1/13  =  304.01  -51.947,  1/14  =  308.71  -93.177,  1715  =  340.82  -  75.327, 

1716  =  259.14-  504.817 

Naiad  ctQ  =  299.36  +  0.70  sin  N  —  6.49  sin  VI  +  0.25  sin  21V 1 
Sg  =  43.36  -  0.51  cos  N  -  4.75  cos  Ml  +  0.09  cos  21V1 
IV  =  254.06  +  1222.8441209 d  -  0.48  sin  N  +  4.40  sin  VI  -  0.27  sin  2N1 
Thalassa  ag  =  299.36  +  0.70  sin  N  —  0.28  sin  N 2 
<>0  =  43.45  -  0.51  cos  N  -  0.21  cos  N2 
W  =  102.06+  1155. 7555612af  —  0.48  sin  N  +  0.19  sin  N2 
Despina  ag  =  299.36  +  0.70sin  N  —  0.09  sin  N 3 
Sg  =  43.45  -  0.51  cos  N  -  0.07  cos  N3 
W  =  306.51  +  1075.7341562 d  -  0.49  sin  N  +  0.06  sin  IV3 
Galatea  ag  =  299.36  +  0.70  sin  N  —  0.07  sin  1V4 
Sg  =  43.43  -  0.51  cos  N  -  0.05  cos  1V4 
IV  =  258.09  +  839.6597686 d  -  0.48  sin  N  +  0.05  sin  N4 
Larissa  ctp  =  299.36  +  0.70  sin  N  —  0.27  sin  N5 
Sg  =  43.41  -  0.51  cos  N  -  0.20  cos  N5 
IV  =  179.41  +  649.0534470(1  -  0.48  sin  N  +  0.19  sin  N5 
Proteus  ag  =  299.27  +  0.70  sin  N  —  0.05  sin  N 6 
<50  =  42.91  -  0.51  cos  N  -  0.04  cos  N6 
IV  =  93.38  +  320.7654228(1  -  0.48  sin  N  +  0.04  sin  N6 
Triton  ag  =  299.36  -  32.35  sin  N1  -  6.28  sin  2N7  -  2.08  sin  3V7 

-  0.74  sin  4N7  -  0.28  sin  5N7  -  0.1 1  sin  6N7 

-  0.07  sin  7N7  -  0.02  sin  8V7  -  0.01  sin  9N7 

So  =  41.17  +  22.55  cos  N1  +  2.10  cos  2N1  +  0.55  cos  3V7 

+  0.16  cos  4N1  +  0.05  cos  5 1V7  +  0.02  cos  6N7 
+  0.01  cos  7N1 
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Table  2  continued 

W  =  296.53  -  61. 2572637c/  +  22.25  sin  N1  +  6.73  sin  2N1 
+2.05  sin  3JV7  +  0.74  sin  4N1  +  0.28  sin  5N1 
+0. 1 1  sin  6V7  +  0.05  sin  7 N1  +  0.02  sin  8 N1 
+0.01  sin  9JV7 

where  N  =  357°.85  +  52°. 3167",  N 1  =  323°. 92  +  62606°  .67.  N2  =  220°.51  +  55064°.27, 

N3  =  354.27  +  46564.577  N4  =  75.31  +  26109.47,  N5  =  35.36  +  14325.47, 

N6  =  142.61  +  2824.67,  N1  =  177.85  +  52.3167 

(a)  Caution:  These  formulae  are  only  precise  to  approximately  150m.  For  higher  precision  an  ephemeris 
should  be  used  as  described  in  Sect.  3  of  the  text 

(b)  The  0°  meridian  is  defined  by  the  sub-Jovian  direction  since  it  is  assumed  surface  features  will  not  last 
long  enough  to  serve  as  a  long  term  reference 

(c)  The  182°  meridian  is  defined  by  the  crater  Cilix 

(d)  The  128°  meridian  is  defined  by  the  crater  Anat 

(e)  The  326°  meridian  is  defined  by  the  crater  Saga 

(f) These  equations  are  correct  for  the  period  of  the  Voyager  encounters.  Because  of  precession  these  may 
change.  Additionally,  orbital  swaps  between  Janus  and  Epimetheus  induce  changes  in  their  mean  spin  rates, 
and  they  are  subject  to  forced  librations 

(g)  The  162°  meridian  is  defined  by  the  crater  Palomides 

(h)  The  5°  meridian  is  defined  by  the  crater  Salih 

(i)  The  299°  meridian  is  defined  by  the  crater  Arete 

(j)  The  63°  meridian  is  defined  by  the  crater  Palinurus 

(k)  The  340°  meridian  is  defined  by  the  crater  Tore 

(l)  The  276°  meridian  is  defined  by  the  crater  Almeric 

Satellites  for  which  no  suitable  data  are  yet  available  have  been  omitted  from  this  table.  Nereid  is  not  included 
in  this  table  because  it  is  not  in  synchronous  rotation 


the  currently  known  orientation  of  the  dynamically  oriented  body  fixed  frame  is  probably 
more  accurately  known  than  that  of  the  established  feature  fixed  frame.  When  such  systems 
and  frames  in  such  systems  are  used,  work  to  derive  the  relationships  between  them  and 
the  recommended  cartographic  coordinate  system  should  be  undertaken  so  that  conversions 
between  the  systems  can  be  accomplished  at  some  known  level  of  accuracy  when  necessary. 
This  will  also  eventually  allow  the  creation  of  final  cartographic  products  in  the  recom¬ 
mended  system.  Users  should  also  be  aware  that  at  high  levels  of  precision  (e.g.  for  the  Moon 
and  probably,  but  not  yet  measured,  for  Mercury)  a  principal  axis  system  is  not  necessarily 
coincident  with  systems  defined  via  principles  of  synchronous  or  resonant  rotation.  Principal 
axis  frames  usually  rely  on  a  specific  gravity  field  model  for  their  definition,  and  may  often 
change  with  improved  gravity  field  determinations — just  as  frames  that  rely  on  fixed  features 
may  often  change  when  a  body  is  first  being  mapped  or  mapped  at  improved  resolution  and 
accuracy. 


3  Coordinate  system  for  the  Moon 

The  recommended  coordinate  system  for  the  Moon  is  the  mean  Earth/polar  axis  (ME)  sys¬ 
tem.  This  is  in  contrast  to  the  principal  axis  (PA)  system,  sometimes  called  the  axis  of  figure 
system.  The  ME  system,  sometimes  called  the  mean  Earth/rotation  axes  system,  is  recom¬ 
mended  because  nearly  all  cartographic  products  of  the  past  and  present  have  been  aligned 
to  it  (Davies  and  Colvin  2000).  The  difference  in  the  coordinates  of  a  point  on  the  surface 
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of  the  Moon  between  these  systems  is  approximately  860  m.  In  past  reports  the  rotation  and 
pole  position  for  the  Moon  have  been  given  for  the  ME  system  using  closed  formulae.  For 
convenience  for  many  users,  those  formulae  are  repeated  here  in  Table  2.  However,  users 
should  note  that  these  are  valid  only  to  the  approximately  150  m  level  of  accuracy,  as  shown 
e.g.,  by  Konopliv  et  al.  (2001,  Figure  3).  For  high  precision  work  involving  e.g.,  spacecraft 
operations,  high-resolution  mapping,  and  gravity  field  determination,  it  is  recommended  that 
a  lunar  ephemeris  be  used  to  obtain  the  libration  angles  for  the  Moon  from  which  the  pole 
position  and  rotation  can  be  derived. 

Specifically,  the  NASA/JPL  planetary  and  lunar  ephemeris  DE  421  (Williams  et  al.  2008; 
Folkner  et  al.  2008,  2009)  is  considered  the  best  currently  available  lunar  ephemeris.  We 
follow  the  recommendations  of  the  Lunar  Reconnaissance  Orbiter  mission  and  the  (NASA) 
Lunar  Geodesy  and  Cartography  Working  Group  (LRO  Project  and  LGCWG  2008)  regard¬ 
ing  its  use.  The  DE  421  ephemeris  may  be  downloaded  in  an  ASCII  version  from  ftp://ssd. 
jpl.nasa.gov/pub/eph/planets/ascii/de421/.  The  development  of  a  new  JPL  lunar  ephemeris 
is  under  consideration  (Williams  2009,  private  communication)  and,  if  it  does  become  avail¬ 
able,  it  might  be  used  for  the  highest  possible  accuracy.  Polynomial  representations  of  the 
(Euler)  lunar  libration  angles  and  their  rates  in  the  PA  system  are  stored  in  the  ephemeris  file. 
These  three  libration  angles  are: 

a)  ip ,  the  angle  along  the  ICRF  equator,  from  the  ICRF  X-axis  to  the  ascending  node  of  the 
lunar  equator; 

b)  9,  the  inclination  of  the  lunar  equator  to  the  ICRF  equator;  and 

c)  \j/,  the  angle  along  the  lunar  equator  from  the  node  to  the  lunar  prime  meridian. 

Coordinates  or  Euler  angles  in  the  ME  system  (vector  M)  can  be  rotated  to  the  PA  system 
(vector  P)  using  the  following  expression  (Williams  et  al.  2008,  p.  10): 

P  =  Rz  (67". 92)  Ry  (78". 56)  Rx  (0".30)  M  (1) 

Conversely,  coordinates  or  Euler  angles  in  the  PA  system  can  be  rotated  into  the  ME  system 
with: 


M  =  Rx  (— 0".30)  Ry  (—78". 56)  Rz  (-67".92)  P  (2) 

where  Rx,  Ry,  and  Rz  are  the  standard  rotation  matrices  for  right-handed  rotations  around 
the  X,  Y,  and  Z  axes  respectively. 

Therefore,  for  a  given  epoch,  the  user  should  obtain  ip,  6,  and  i/r  from  the  ephemeris  file 
and  store  them  as  the  vector  P,  apply  the  transformation  in  Eq.  2,  and  extract  the  angles,  now 
in  the  ME  system,  from  the  vector  M.  These  angles  can  then  be  converted  with: 

ao  =  <p  —  90° 

<50  =  90°  -  9 
W  =  i/r 

giving  the  lunar  rotation  angles  in  the  standard  ao,  and  W  formulation  (of  Table  2)  and  in 
the  ME  system. 

Alternatively,  if  the  user  has  coordinates  for  a  point  in  ICRF  coordinates  (vector  I)  that 
they  wish  to  convert  to  ME  coordinates,  for  a  given  epoch  the  user  should  obtain  ip,  6,  and  i/r 
from  the  ephemeris  file,  and  then  do  the  conversion: 

M  =  Rx  (— 0".30)  Ry  (—78". 56)  Rz  (-67".92)  Rz  (t/0  Rx  (9)  Rz  (ip)  I  (3) 
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with  M  now  being  the  coordinates  of  the  point  in  the  ME  system.  The  user  should  note  that 
the  numerical  values  for  the  rotations  in  Eqs.  1,  2,  and  3  are  specific  to  DE  421  and  are 
different  for  past  and  future  ephemerides. 

Note  that  the  NASA/JPL  Navigation  and  Ancillary  Information  Facility  (NAIF)  provides 
software  and  files  to  facilitate  the  above  transformations.  This  includes  a  Planetary  Con¬ 
stants  Kernel  (PCK)  made  using  the  lunar  libration  information  extracted  from  the  DE  42 1 
ephemeris,  and  a  special  lunar  frames  kernel  (FK)  providing  the  specifications  and  data 
needed  to  construct  the  PA  to  ME  system  transformation.  A  new  version  of  the  PCK  will  also 
be  provided  when  a  new  IPL  ephemeris  is  released.  See  http://naif.jpl.nasa.gov  or  (for  DE 
421)  http://naif.jpl.nasa.gov/naif/lunar_kernels.txt  for  further  information.  Although  written 
before  DE  421  became  available,  Roncoli  (2005)  also  provides  useful  information  on  lunar 
constants  and  coordinates. 


4  Rotation  elements  for  planets  and  satellites 

The  orientation  model  for  Mercury  has  been  updated  using  improved  pole  position  and 
short  period  longitude  libration  information  from  Margot  (2009).  Margot  also  recommends 
(in  essence)  returning  the  longitude  origin  of  Mercury  to  a  dynamically  defined  one,  ori¬ 
ented  along  the  long  principal  axis  of  Mercury,  as  it  was  in  the  past.  However,  following  the 
practice  of  this  Working  Group  since  its  original  report  (Davies  et  al.  1980)  and  the  established 
principle  that  “where  possible,  however,  the  cartographic  position  of  the  prime  meridian  is 
defined  by  a  suitable  observable  feature,”  we  continue  to  use  the  crater  Hun  Kal  (which  means 
“twenty”  in  the  Mayan  language)  to  define  the  20°  meridian.  Therefore,  we  base  the  value  of 
Wo  on  that  used  previously  (Wo  =  329°. 548,  Robinson  et  al.  1999;  Seidelmann  et  al.  2002, 
2005,  2007),  although  corrected  for  the  total  of  the  libration  terms  at  J2000.0  (y./2000.o)-  This 
results  in  a  new  recommended  value  of: 

Wo  =  329°. 548  -  x,2ooo.o  =  329°. 5469 

We  encourage  and  look  forward  to  a  new  determination  of  this  value  from  new  MESSENGER 
or  a  combination  of  MESSENGER  and  Mariner  10  mapping. 

The  rotation  rate  of  Saturn,  that  is  given  in  Table  1  is  based  on  Voyager  observations  of 
kilometer  wavelength  radio  signals.  Recent  Cassini  observations  (Giampieri  et  al.  2006)  of  a 
signal  in  Saturn’s  magnetic  field  gives  a  period  of  about  lOh  and  47  min,  about  8  min  longer 
than  the  previously  determined  period.  At  this  time,  it  is  still  uncertain  whether  this  is  the  true 
rotation  rate  or  what  physical  mechanism  is  causing  the  different  signals  (Stevenson  2006). 
See  Kurth  et  al.  (2007),  Gurnett  et  al.  (2007),  Anderson  and  Schubert  (2007),  and  Russell 
and  Dougherty  (2010)  for  additional  discussion.  Russell  reports  (2010,  private  communica¬ 
tion)  though  that  the  Russell  and  Dougherty  (2010)  results  indicating  a  specific  period  are 
in  doubt,  and  it  appears  as  if  the  rotation  period  cannot  be  obtained  reliably  from  Cassini 
magnetic  field  observations  alone.  He  and  others  involved  with  the  Cassini  Project  (Spilker 
2010,  private  communication)  are  looking  into  whether  some  general  or  consensus  solution 
can  be  achieved.  Hence,  the  rotation  rate  of  Saturn  has  not  been  changed,  while  more  results 
from  the  Cassini  mission  are  anticipated. 

The  rotation  rates  of  Uranus  and  Neptune  were  determined  from  the  Voyager  mission  in 
1986  and  1989.  The  uncertainty  of  those  rotation  rates  are  such  that  the  uncertainty  of  the 
actual  rotation  position  in  the  present  day  is  more  than  a  complete  rotation  in  each  case. 

A  new  model  for  the  pole  position  and  rotation  of  Mars  has  been  proposed  by  Konopliv 
et  al.  (2006)  based  on  the  most  recent  spacecraft  data.  At  this  time,  following  the  advice  of 
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the  NASA  Mars  Geodesy  and  Cartography  Working  Group  and  its  chair  (Duxbury  2006; 
Duxbury  2009,  private  communication),  the  use  of  this  new  model  is  not  recommended  for 
cartographic  purposes.  This  is  for  a  number  of  reasons  including  that  for  the  immediate  future 
the  new  model  would  have  little  if  any  effect  on  cartographic  products,  and  also  that  it  is 
expected  to  be  significantly  changed  in  the  next  few  years  as  new  data  become  available. 
However,  users  with  high  accuracy  requirements,  such  as  Mars  gravity  field  determination, 
may  wish  to  consider  using  it. 


5  Rotational  elements  for  dwarf  planets,  minor  planets,  their  satellites,  and  comets 

For  planets  and  satellites,  the  IAU  definition  of  north  pole  is  the  pole  that  lies  above  the 
invariant  plane  of  the  Solar  System,  and  the  rotation  can  be  either  direct  or  retrograde.  For 
dwarf  planets,  minor  planets,  their  satellites,  and  comets,  given  substantial  indirect  evidence 
for  large  precession  of  the  rotational  poles  of  some  comets,  this  first  definition  needs  to  be 
rethought.  In  particular,  situations  exist  in  which  the  pole  that  is  clearly  “north”  in  the  IAU 
sense  precesses  over  several  decades  to  become  clearly  “south”  in  the  IAU  sense.  Comet 
2P/Encke  is  a  prime  example  of  a  comet  for  which  very  large  precession  has  been  inferred. 

There  is  also  clear  evidence  for  excited  state  rotation  for  comet  1  P/Halley  and  minor 
planet  (4179)  Toutatis.  In  this  case,  the  angular  momentum  vector  moves  around  on  the  sur¬ 
face  of  the  body.  The  rotational  spin  vector  describes  substantial  excursions  from  the  angular 
momentum  vector  during  the  course  of  the  7-day  periodicity  that  is  seen  in  the  light  curve. 
We  can,  therefore,  anticipate  cases  in  which  the  rotational  spin  pole  moves  back  and  forth 
between  north  and  south  on  a  time  scale  of  days.  Thus,  there  is  the  issue  of  needing  to  change 
our  definition  of  the  rotational  pole. 

The  choice  of  a  rotational  pole  for  a  body  in  simple  rotation  with  slow  precession  is 
straightforward.  One  can  choose  the  pole  that  follows  either  the  right-hand  rule  or  the  left- 
hand  rule,  and  the  right-hand  rule  is  chosen  here.  This  would  be  the  “positive”  pole  to  avoid 
confusion  with  the  north-south  terminology.  Ideally  one  would  like  to  choose  a  pole  for 
excited  state  rotation  that  reduces  to  this  definition  as  the  rotational  energy  relaxes  to  the 
ground  state.  For  SAM  (short-axis  mode)  rotational  states,  it  is  possible  to  define  a  body- 
fixed  axis  that  circulates  in  a  generally  complex  pattern  about  the  angular  momentum  vector 
and  this  approaches  the  simple  right-hand  rule  definition  as  the  rotational  energy  relaxes  to 
the  ground  state  of  simple  rotation.  Presumably  the  appropriate  body-fixed  pole  is  the  axis  of 
maximum  moment  of  inertia.  However,  the  definition  for  a  body  in  a  LAM  (long-axis  mode) 
rotational  state  is  not  so  obvious,  because  there  is  then  complete  rotation  about  the  long  axis 
of  the  body  as  well  as  rotation  about  a  short  axis.  In  this  case,  the  pole  should  be  taken  as  the 
minimum  moment  of  inertia  (the  long  axis  of  an  ellipsoid)  according  to  the  right-hand  rule. 

In  practice,  the  initial  encounter  with  a  small  irregular  body  may  not  provide  enough  infor¬ 
mation  to  determine  the  shape,  moments  of  inertia,  and  rotational  dynamics  with  sufficient 
accuracy  that  rotational  parameters  based  on  them  will  stand  the  test  of  time.  In  such  cases, 
the  recommended  approach  to  defining  rotational  parameters  and  coordinate  axes  should 
be  based  on  the  same  general  principles  that  apply  to  planets  and  satellites.  If  possible,  the 
initial  definition  of  a  body-fixed  coordinate  system  should  be  based  on  a  shape  model  and 
estimate  of  the  moments  of  inertia,  with  the  polar  axis  chosen  as  described  above.  If  there  is 
insufficient  information  to  determine  the  moments  of  inertia,  it  may  be  necessary  to  define 
a  coordinate  system  based  on  the  instantaneous  axis  of  rotation  at  the  time  of  encounter. 
The  choice  of  prime  meridian  is  potentially  arbitrary,  but  there  is  precedent  (e.g.  with  (433) 
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Eros)  for  choosing  it  so  it  aligns  with  the  longest  axis  (or  minimum  moment  of  inertia,  if  this 
can  be  estimated). 

The  orientation  of  the  body  at  subsequent  encounters  is  likely  to  differ  from  the  predic¬ 
tions  of  a  model  based  on  initial  encounter  data.  Such  departures  may  be  gross,  as  a  result 
of  rapid  precession  in  which  the  spin  axis  varies  relative  to  the  body,  or  more  subtle,  simply 
because  of  the  limits  of  accuracy  of  the  initial  observation.  In  either  case,  it  should  be  borne 
in  mind  that  (as  for  planets  and  satellites)  the  main  purpose  of  defining  a  body-fixed  coor¬ 
dinate  system  is  to  facilitate  mapping  of  surface  features.  It  is,  therefore,  desirable  to  relate 
the  axes  of  the  initial  system  to  identifiable  surface  features.  When  new  observations  become 
available,  the  axes  should  in  most  cases  be  left  unchanged  with  respect  to  the  surface  features 
and  the  rotational  model  amended  to  model  the  inertial-space  orientation  of  the  axes  more 
accurately.  In  contrast  to  the  case  of  planets  and  satellites,  for  which  the  rotation  axis  is  often 
determined  more  accurately  than  the  rotation  rate,  for  irregular  bodies  both  axis  and  rate  may 
vary  over  time  or  be  poorly  determined.  It  follows  that  two  or  three  landmark  features  will 
be  required  to  determine  the  body-fixed  orientation  of  the  coordinate  axes,  rather  than  the 
single  landmark  that  suffices  to  define  the  prime  meridian  of  a  regularly  rotating  body. 

As  specified  in  Sect.  6,  for  planets  and  satellites,  longitude  should  increase  monotoni- 
cally  for  an  observer  fixed  in  inertial  space.  For  dwarf  planets,  minor  planets,  their  satellites, 
and  comets  however,  with  the  above  rule  for  poles,  this  definition  corresponds  always  to  a 
left-hand  rule  for  increasing  longitude,  since  the  concept  of  retrograde  rotation  is  no  longer 
relevant.  Therefore,  for  dwarf  planets,  minor  planets,  their  satellites,  and  comets,  to  be  consis¬ 
tent  with  the  above  pole  definition,  increasing  longitude  should  always  follow  the  right-hand 
rule.  This  definition  is  consistent  with  the  sense  of  increasing  longitude  used  for  Eros  by 
Miller  et  al.  (2002),  but  is  inconsistent  with  the  sense  of  increasing  longitude  used  for  Eros 
by  Thomas  et  al.  (2002). 

For  each  such  body,  the  positive  pole  of  rotation  is  selected  as  the  maximum  or  minimum 
moment  of  inertia  according  to  whether  there  is  short  or  long  axis  rotational  state  and  accord¬ 
ing  to  the  right-hand  rule.  So  the  positive  pole  is  specified  by  the  value  of  its  right  ascension 
ao  and  declination  <5o-  With  the  pole  so  specified,  the  two  intersection  points  of  the  body’s 
equator  and  the  ICRF  equator  are  ao  ±  90°.  We  choose  one  of  these,  ao  +  90°,  and  define 
it  as  the  node  Q.  Suppose  the  prime  meridian  has  been  chosen  so  that  it  crosses  the  body’s 
equator  at  the  point  B.  We  then  specify  the  location  of  the  prime  meridian  by  providing 
a  value  for  W,  the  angle  measured  along  the  body’s  equator  between  the  node  Q  and  the 
point  B  in  a  right-hand  system  with  respect  to  the  body’s  positive  pole  (see  Fig.  2).  The  right 
ascension  of  the  point  Q  is  90°  +  ao  and  the  inclination  of  the  body’s  equator  to  the  celestial 
equator  is  90°  —  8q.  As  long  as  the  planet,  and  hence  its  prime  meridian,  rotates  uniformly, 
W  varies  linearly  with  time  according  to  the  right-hand  rule.  In  addition,  ao, ,  So,  and  W  may 
vary  with  time  due  to  a  precession  of  the  axis  of  rotation  of  the  body.  It  should  be  noted  that 
for  bodies  whose  spin  precesses  rapidly  with  large  amplitude,  this  simple  formulation  of  the 
body  orientation  in  terms  of  pole  orientation  and  spin  angle  may  be  insufficient.  In  such  cases 
it  may  be  necessary  to  develop  expressions  for  the  body  orientation  in  terms  of  a  full  set  of 
time-varying  Euler  angles.  At  the  present  time,  however,  the  formulation  given  here  suffices 
to  represent  the  rotation  of  those  bodies  for  which  data  are  available. 

The  angle  W  specifies  the  ephemeris  position  of  the  prime  meridian,  and  for  dwarf  planets, 
minor  planets,  their  satellites,  and  comets  without  any  accurately  observable  fixed  surface 
features,  the  adopted  expression  for  W  defines  the  prime  meridian.  Where  possible,  however, 
the  cartographic  position  of  the  prime  meridian  is  defined  by  a  suitable  observable  feature, 
and  so  the  constants  in  the  expression  W  =  Wq  +  W  d,  where  d  is  the  interval  in  days  from 
the  standard  epoch,  are  chosen  so  that  the  ephemeris  position  follows  the  motion  of  the 
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Fig.  2  Reference  system  used  to 
define  orientation  of  dwarf 
planets,  minor  planets,  their 
satellites,  and  comets 


cartographic  position  as  closely  as  possible;  in  these  cases  the  expression  for  W  may  require 
emendation  in  the  future.  Table  3  gives  the  recommended  rotation  values  for  the  direction  of 
the  positive  pole  of  rotation  and  the  prime  meridian  of  selected  dwarf  planets,  minor  planets, 
their  satellites,  and  comets.  Values  are  given  for  objects  that  have  been  imaged  by  spacecraft, 
radar,  or  high  resolution  Earth-based  imaging  systems  with  sufficient  resolution  to  estab¬ 
lish  accurate  pole  orientation  and  rotation  rates.  Values  are  not  given  for  objects  where  the 
observations  are  limited  to  photometric  light  curves. 

When  new  higher  accuracy  mapping  is  done,  the  longitudes  and  possibly  latitudes  of  the 
fixed  feature  or  features  should  be  maintained  and  a  new  value  for  Wo  (and,  as  necessary, 
the  other  rotational  parameters)  derived,  published  in  peer-reviewed  literature,  and  reported 
to  this  Working  Group  for  possible  adoption.  For  bodies  where  they  are  in  use,  longitude 
defining  features  are  noted  in  the  footnotes  to  Table  3.  As  already  described  in  Sect.  2  above, 
once  one  or  more  features  has  or  have  been  chosen  with  defined  cartographic  coordinates, 
they  should  not  be  changed  except  in  unusual  circumstances  or  in  cases  where  the  change 
amounts  to  refining  the  landmark  position  within  the  precision  to  which  it  was  originally 
measured. 

A  dwarf  planet,  minor  planet,  one  of  their  satellites,  or  a  comet  will  typically  be  included 
in  Table  3  only  if  it  meets  fundamental  publication,  data  quality,  and  applicability  criteria. 
Estimated  values  for  ao>  <5o.  and  Wo  should  appear  in  a  publication  that  appears  in  a  refereed 
journal  (or  a  revision  that  is  in  preparation).  The  analysis  to  determine  these  values  will  have 
been  derived  from  data  of  sufficient  fidelity  and  quality  to  assure  an  accurate  estimate,  and 
some  portion  of  that  data  will  have  been  acquired  via  direct  methods  (e.g.,  direct  imaging 
from  a  spacecraft,  a  space  telescope,  or  an  adaptive  optics  system)  although  possibly  used  in 
combination  with  data  from  other  methods,  including  indirect  ones  (e.g.  photometry,  multi¬ 
chord  stellar  occultation).  Lastly,  a  cartographic  need  must  exist  that  justifies  the  definition 
of  a  prime  meridian  and  pole  for  this  body. 

As  described  above  for  planets  and  satellites,  these  recommendations  are  not  intended  to 
imply  that  other  coordinate  systems  with  different  rotational  elements  should  not  be  used  for 
minor  planets  and  comets  for  other  than  cartographic  purposes.  See  Sect.  2  for  a  discussion  of 
possible  options,  such  as  using  a  principal  axis  system  for  other  than  cartographic  purposes. 
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Table  3  Recommended  rotation  values  for  the  direction  of  the  positive  pole  of  rotation  and  the  prime  meridian 
of  selected  dwarf  planets,  minor  planets,  their  satellites,  and  comets 

d  is  the  interval  in  days  from  the  standard  epoch,  i.e.  J2000.0  =  JD 

2451545.0,  i.e.  2000  January  1  12 hours  TDB.  g?o,  <5o>  and  Ware  as  defined 
in  the  text 

(1)  Ceres 

<*0 

=  291°  ±5° 

<5o 

=  59°  ±  5° 

w 

=  170°  .90  +  952°.1532rf 

(a) 

(2)  Pallas 

01 0 

=  33° 

So 

=  -3° 

w 

=  38°  +  1 105°.  8036rf 

(b) 

(4)  Vesta 

“0 

=  305°. 8  ±  3°.l 

80 

=  41°.4d=  1°.5 

w 

=  292°+  1617°.332776c( 

(c) 

(21)  Lutetia 

Oi  0 

=  52°  ±  5° 

So 

=  12°  ±  5° 

w 

=  94°  +  1057°.7515rf 

(d) 

(243)Ida 

“0 

=  168°. 76 

=  -2°. 88 

w 

=  265°  .95  +  1864°. 6280070(7 

(e) 

(433)  Eros 

ao 

=  11°.35±0°.02 

So 

=  17°. 22  ±  0°.02 

w 

=  326°  .07+  1639°.38864745r/ 

(f) 

(511)  Davida 

ao 

=  297° 

80 

=  5° 

w 

=  268°. 1  +  1684°  .4193549,7 

(g)  (h) 

(95 1 )  Gaspra 

»0 

=  9°  .47 

So 

=  26°. 70 

w 

=  83°. 67  +  1226°. 91 14850,7 

(i) 

(2867)  Steins 

ao 

=  90° 

80 

=  -62° 

w 

=  93°. 94  +  1428°. 852332,7 

(j) 

(25143)  Itokawa 

“0 

=  90°. 53 

So 

=  -66°.  30 

w 

=  000°  +  712°. 143d 

(k) 

(134340)  Pluto 

ao 

=  132°. 993 

So 

=  -6°. 163 

w 

=  237°. 305  +  56°. 3625225,7 

(1) 

(134340)  Pluto  :  I  Charon 

“0 

=  132°. 993 

80 

=  -6°. 163 

w 

=  57°. 305  +  56°. 3625225,7 

(m) 

9P/Tempel  1 

ao 

=  294° 

So 

=  73° 

w 

=  252°. 63  +  2 12°  .064,7 

(n) 
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Table  3  continued 

d  is  the  interval  in  days  from  the  standard  epoch,  i.e.  J2000.0  =  JD 
2451545.0,  i.e.  2000  January  1  12 hours  TDB.  c^q,  <5q,  and  Ware  as  defined 
in  the  text 


19P/Borrelly  a0  =  2 18°. 5  ±  3° 

So  = -12°.5  ±  3° 

W  =  000°  +  390°  ,0d  (k) 

(a)  The  0°  meridian  is  defined  by  the  unnamed  bright  spot  shown  in  Fig.  1  in  Thomas  et  al.  (2005)  and  Figures  5 
and  6  at  the  0°  meridian  in  Li  et  al.  (2006a) 

(b)  The  0°  meridian  is  defined  by  the  direction  (positive  x)  of  the  long  axis  of  the  Carry  et  al.  (2010a) 
shape  model 

(c)  The  0°  meridian  is  defined  by  a  feature  not  yet  formally  named,  but  referred  to  as  Olbers  Regio  by 
Thomas  et  al.  (1997) 

(d)  The  0°  meridian  has  (so  far)  been  arbitrarily  defined  based  on  light  curve  information 

(e)  The  0°  meridian  is  defined  by  the  crater  Afon 

(f)  The  0°  meridian  is  defined  by  an  unnamed  crater 

(g)  The  0°  meridian  is  defined  by  the  direction  of  the  long  axis  that  points  toward  the  Earth  on  2002  December 

"  27  7.83  UT  (Conrad  et  al.  2007) 

(h)  Values  have  been  revised  from  those  which  appear  in  Conrad  et  al.  (2007 )  to  the  values  given  above,  which 
appear  in  a  publication  by  the  same  authors  (in  preparation) 

(i)  The  0°  meridian  is  defined  by  the  crater  Charax 

(j)  The  0°  meridian  is  defined  by  a  feature  not  yet  formally  named,  but  referred  to  as  the  crater  Spinel  by 
Jorda  et  al.  (2010) 

(k)  Since  only  rotation  rate  information  is  available,  the  0°  meridian  is  currently  arbitrarily  defined  with 
W0  =  0° 

(l)  The  0°  meridian  is  defined  as  the  mean  sub-Charon  meridian 

(m)  The  0°  meridian  is  defined  as  the  mean  sub-Pluto  meridian 

(n)  The  0°  meridian  is  defined  by  a  350  m  diameter  unnamed  circular  feature  near  the  Deep  Impactor  impact 
site  (Thomas  et  al.  2007).  We  expect  that  the  rotational  elements  will  be  changed  in  our  next  report  to  show 
a  time-variable  rotation  rate 


6  Definition  of  cartographic  coordinate  systems  for  planets  and  satellites 

In  mathematical  and  geodetic  terminology,  the  terms  ‘latitude’  and  ‘longitude’  refer  to  a 
right-hand  spherical  coordinate  system  in  which  latitude  is  defined  as  the  angle  between  a 
vector  passing  through  the  origin  of  the  spherical  coordinate  system  and  the  equator,  and 
longitude  is  the  angle  between  the  vector  and  the  plane  of  the  prime  meridian  measured  in 
an  eastern  direction.  This  coordinate  system,  together  with  Cartesian  coordinates,  is  used  in 
most  planetary  computations,  and  is  sometimes  called  the  planetocentric  coordinate  system. 
In  this  system,  longitudes  are  always  positive  toward  the  east.  The  origin  is  the  center  of  mass. 

Because  of  astronomical  tradition,  planetographic  coordinates  (those  commonly  used  on 
maps)  may  or  may  not  be  identical  with  traditional  spherical  coordinates.  Planetographic 
coordinates  are  defined  by  guiding  principles  contained  in  a  resolution  passed  at  the  four¬ 
teenth  General  Assembly  of  the  IAU  in  1970.  These  guiding  principles  state  that: 

(1)  The  rotational  pole  of  a  planet  or  satellite  which  lies  on  the  north  side  of  the  invariable 
plane  will  be  called  north,  and  northern  latitudes  will  be  designated  as  positive. 

(2)  The  planetographic  longitude  of  the  central  meridian,  as  observed  from  a  direction  fixed 
with  respect  to  an  inertial  system,  will  increase  with  time.  The  range  of  longitudes  shall 
extend  from  0°  to  360°. 
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Thus,  west  longitudes  (i.e.,  longitudes  measured  positively  to  the  west)  will  be  used  when 
the  rotation  is  direct,  i.e.  the  sign  of  the  second  term  in  the  expression  for  W  is  positive.  East 
longitudes  (i.e.,  longitudes  measured  positively  to  the  east)  will  be  used  when  the  rotation  is 
retrograde,  i.e.  the  sign  of  the  second  term  in  the  expression  for  W  is  negative.  The  origin  is 
the  center  of  mass.  Also  because  of  tradition,  the  Earth,  Sun,  and  Moon  do  not  conform  with 
this  definition.  Their  rotations  are  direct  and  longitudes  run  both  east  and  west  180°,  or  east 
360°. 

For  planets  and  satellites,  latitude  is  measured  north  and  south  of  the  equator;  north  lat¬ 
itudes  are  designated  as  positive.  The  planetographic  latitude  of  a  point  on  the  reference 
surface  is  the  angle  between  the  equatorial  plane  and  the  normal  to  the  reference  surface 
at  the  point.  In  the  planetographic  system,  the  position  of  a  point  (P)  not  on  the  reference 
surface  is  specified  by  the  planetographic  latitude  of  the  point  (P')  on  the  reference  surface 
at  which  the  normal  passes  through  P  and  by  the  height  (h)  of  P  above  P' . 

The  topographic  reference  surface  of  Mars  is  that  specified  in  the  final  MOLA  Mission 
Experiment  Gridded  Data  Record  (MEGDR)  Products  (Smith  et  al.  2003).  In  particular, 
the  128  pixels/0  resolution,  radius  and  topographic  surfaces  are  recommended,  although  the 
lower  resolution  versions  may  be  used  where  appropriate  and  documented,  and  for  the  areas 
poleward  of  ±88°  latitude. 

For  Mercury,  the  use  of  a  planetocentric,  east-positive  (right-handed)  system  was  adopted 
by  the  MESSENGER  project  more  than  9  years  ago  to  facilitate  geodetic  analysis,  particu¬ 
larly  topography  and  gravity,  as  well  as  all  cartography.  The  Mariner  10  mission  used  the 
IAU  standard  system.  There  are  standard  transformations  between  the  two  coordinate  sets. 
For  the  Mars  Global  Surveyor  mission,  an  areocentric,  east-positive  system  was  used  despite 
years  of  Mariner  4,  6,  7,  and  9  and  Viking  data  mapped  with  the  IAU  standard  system. 

The  reference  surfaces  for  some  planets  (such  as  Earth  and  Mars)  are  ellipsoids  of  revo¬ 
lution  for  which  the  radius  at  the  equator  (A)  is  larger  than  the  polar  semi-axis  (C). 

Calculations  of  the  hydrostatic  shapes  of  some  of  the  satellites  (Io,  Mimas,  Enceladus,  and 
Miranda)  indicate  that  their  reference  surfaces  should  be  triaxial  ellipsoids.  Triaxial  ellip¬ 
soids  would  render  many  computations  more  complicated,  especially  those  related  to  map 
projections.  It  would  be  difficult  to  generalize  many  projections  so  as  to  retain  their  elegant 
and  popular  properties  and  there  is  a  lack  of  agreement  on  basic  matters  such  as  the  appro¬ 
priate  definitions  of  latitude  and  longitude.  For  these  reasons  spherical  reference  surfaces  are 
frequently  used  in  mapping  programs. 

Many  small  bodies  of  the  Solar  System  (satellites,  minor  planets,  and  comet  nuclei)  have 
very  irregular  shapes.  Sometimes  spherical  reference  surfaces  are  used  for  computational 
convenience,  but  this  approach  does  not  preserve  the  area  or  shape  characteristics  of  com¬ 
mon  map  projections.  Orthographic  projections  often  are  adopted  for  cartographic  portrayal 
as  these  preserve  the  irregular  appearance  of  the  body  without  artificial  distortion.  A  more 
detailed  discussion  of  cartographic  coordinate  systems  for  small  bodies  is  given  in  Sect.  7  of 
this  report. 

Table  4  gives  the  size  and  shape  parameters  for  the  planets.  The  Sun  is  included  for  com¬ 
parison  purposes.  Average  (AVG),  north  (N),  and  south  (S)  polar  radii  are  given  for  Mars. 
For  the  purpose  of  adopting  a  best-fitting  ellipsoid  for  Mars,  the  average  polar  radius  should 
be  used — the  other  values  are  for  comparison  only,  e.g.  to  illustrate  the  large  dichotomy  in 
shape  between  the  northern  and  southern  hemispheres  of  Mars.  In  applications  where  these 
differences  may  cause  problems,  the  earlier  recommended  topographic  shape  model  for  Mars 
should  probably  be  used  as  a  reference  surface.  The  mean  radii  shown  in  Tables  4,  5,  and  6  in 
general  (the  Earth  is  an  exception)  are  from  the  original  authors  and  have  not  been  computed 
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from  the  other  radii  by  the  Working  Group  on  the  assumption  that  at  least  some  of  them  are 
independently  computed. 

Table  5  gives  the  size  and  shape  of  satellites  where  known.  Only  brightnesses  are  known 
for  many  of  the  newly  discovered  satellites.  Poles  and  rotation  rates  are  also  not  yet  known  for 
the  new  discoveries,  so  those  satellites  are  not  listed.  A  mean  radius  and  best  fitting  triaxial 
radii  are  given  for  Titan.  The  triaxial  radii  are  for  comparison  only,  as  the  differences  are 
small  enough  that  the  mean  radius  should  be  used  for  most  cartographic  purposes. 

The  values  of  the  radii  and  axes  in  Tables  4  and  5  are  derived  by  various  methods  and  do 
not  always  refer  to  common  definitions.  Some  use  star  or  spacecraft  occultation  measure¬ 
ments,  some  use  limb  fitting,  others  use  altimetry  measurements  from  orbiting  spacecraft, 
and  some  use  control  network  computations.  For  the  Earth,  the  spheroid  refers  to  mean  sea 
level,  clearly  a  very  different  definition  from  other  bodies  in  the  Solar  System. 

The  uncertainties  in  the  values  for  the  radii  and  axes  in  Tables  4  and  5  are  generally  those 
of  the  authors,  and,  as  such,  frequently  have  different  meanings.  Sometimes  they  are  standard 
errors  of  a  particular  data  set,  sometimes  simply  an  estimate  or  expression  of  confidence. 
The  radii  and  axes  of  the  large  gaseous  planets,  Jupiter,  Saturn,  Uranus,  and  Neptune  in 
Table  4  refer  to  a  one-bar-pressure  surface.  The  radii  given  in  the  tables  are  not  necessarily 
the  appropriate  values  to  be  used  in  dynamical  studies;  the  radius  actually  used  to  derive  a 
value  of  J2  (for  example)  should  always  be  used  in  conjunction  with  it.  In  Table  5,  ellipsoidal 
fit  axes  of  objects  less  than  200  km  in  radius  are  for  convenient  comparison  and  their  use  for 
any  modeling  can  only  be  approximate. 


7  Cartographic  coordinates  for  dwarf  planets,  minor  planets,  their  satellites, 
and  comets 

For  large  bodies,  a  spherical  or  ellipsoidal  model  shape  has  traditionally  been  defined  for 
mapping,  as  in  our  past  reports.  For  irregularly  shaped  bodies  the  ellipsoid  is  obviously 
useless,  except  perhaps  for  dynamical  studies.  For  very  irregular  bodies,  the  concept  of  a  ref¬ 
erence  ellipsoid  ceases  to  be  useful  for  most  purposes.  For  these  bodies,  topographic  shapes 
are  usually  represented  by  a  grid  of  radii  to  the  surface  as  a  function  of  planetocentric  latitude 
and  longitude  (when  possible,  or  also  by  a  set  of  vertices  and  polygons). 

Another  problem  with  smaller  bodies  is  that  two  coordinates  (i.e.  spherical  angular  mea¬ 
sures)  may  not  uniquely  identify  a  point  on  the  surface  of  the  body.  In  other  words  it  is 
possible  to  have  a  line  from  the  center  of  the  object  intersect  the  surface  more  than  once. 
This  can  happen  on  large  and  even  mostly  ellipsoidal  objects  such  as  the  Earth,  because  of 
such  features  as  overhanging  cliffs  and  natural  bridges  and  arches.  However,  on  large  bodies 
these  features  are  relatively  very  small  and  often  ignored  at  the  scale  of  most  topographic 
maps.  For  small  bodies  they  may  be  fairly  large  relative  to  the  size  of  the  body.  Example 
cases  are  on  Eros  (at  a  small  patch  west  of  Psyche),  and  certainly  on  Kleopatra  (Ostro  et  al. 
2000),  possibly  on  Toutatis  near  its  ‘neck,’  and  perhaps  near  the  south  pole  of  Ida,  some 
radii  may  intersect  the  surface  more  than  once.  Even  on  small  bodies  this  problem  is  usually 
restricted  to  small  areas,  but  it  still  may  make  a  planetocentric  coordinate  system  difficult 
to  use.  Cartographers  always  have  ad  hoc  tricks  for  a  specific  map,  such  as  interpolating 
across  the  problem  area  from  areas  which  are  uniquely  defined,  or  by  showing  overlapping 
contours.  A  Cartesian  or  other  coordinate  geometry  may  be  preferable  for  arbitrarily  complex 
shapes,  such  as  a  toroidal  comet  nucleus,  where  an  active  region  has  eaten  its  way  through 
the  nucleus.  Such  coordinate  geometries  may  also  be  useful  for  irregular  bodies  imaged  only 
on  one  side,  such  as  for  19P/Borrelly  and  81P/Wild  2. 


Springer 


Report  of  the  IAU  Working  Group 


123 


.a 

£ 


52  u 

S  -S 

2  & 


sq 

■'t 


CM  +1 
»o  in 


iq  O 
co  co 


c  ^ 
c  m 

03  co 


P  S  -H 
-H  -H  g 

tq  GO  O 
Os  — I  — I 


S.U  3 

>  W 


"H  -H 


a 


o 

-H 


C_G  =  ~  oo  so 

so  5  £  r-  os 

Os  c3  c3  CO  CO 

sO  CO  GO  SO  CO 


-H  -H  -H  -H 

*  3  ^  - 


^ 

-H  -H  -H 

04  OO  Os 
os  so  in 
'sT  04  in 
O  ■  “ 


SO  04  04 


-H  -H  -H  -H 


tD  Z 


^  Springer 


Table  5  Size  and  shape  parameters  of  the  satellites 


124 


B.  A.  Archinal  et  al. 


cn 


in 


~o  <u 


in  »n  cn 
<N  O  O 


CN 

cn 


cn  so 

d  d 


3 


6 


<U 

00 


~  r"  cn 

a  in  H  h  to 


r-; 

in 

oo 


<u 

£ 


c3 

00 


u 

E 


c3 

00 


s  -c  e 

bJ)  O  on 
C  03  2 

,2  =  ^ 

- - 3  O^  jrt 
<.  <U  V- 


S  '§  ~ 

g.  a  g 

■S 

oo  <0  2 


E  E 

03  c3 

00  00 


£  © 
c3  cn 
cn 


-H 


-sf  SO 
Os  <N 


cn 


£  E 


in 


-H 


cn 

-H  -H 

in  cn 


O  — I  —I 

-H  -H  -H 

oo  cn  cn 

d  — '  d 

'  “  cn  ’—i 


so  CN  00  00 


os  oo  in  so 


cn  cn  in  oo 


o  2 
S  73 
S  < 
>  > 
X  X 


-C 

H 

> 


U 

> 


B  1 

X  > 


o 

-H 

o 


■—  2S 


ffl  < 


£ 

§ 

U 


>  X  X 


tS 


4!)  Springer 


VIII  Pasiphae 

IX  Sinope 


Table  5  continued 


Report  of  the  IAU  Working  Group 


125 


D-  S 
<L>  M 
TD  w 


s 

22 


£  53 
S  13 


D-  -O 


o  5 

Dh  w 


Jo 

?3  'c' 

u  j2 
c  -c  C, 

«  O  W 
O.  C3  2 

5  ?  2 


6 

22 


oo  oo 
o  o 
-H  -H 

<N  ^ 

d  d  oi  ro  in 


o 

-H 


d 

3 

Oh 


-H  -H  -H  -H 

t"~  o  — <  <n 


ra 

22 

Os 

o 

d 

<N 

O 

-H 

r-; 

d 

-H 

-H 

-H 

-H 

in 

cr 

<u 

TD 

c3 

in 

d 

d 

Os 

ro 

o 

■"t 

OJ  1-H  00  o 

--  d  -2  oi 

-H  -H  -H  -H 

OO  O  On 

d  d  d  d 

Ol  O  <n  SO 


ov  in  oo 

d  oi  d  d 

-H  -H  -H  -H 


<  M  a 

^  HH  K  pj 

>  >  >  „ 

X  X  X  X 


-H 

SO 


O 

-H 


-H  -H 


d 

-H 


-H 

so 


i 

i 


-'t 

d 

-H 


-H  -H 


oi  — < 

d  -2 
-H  -H 

co 

d  00 
in  <N 


-H  -H 

O)  ^t 


N  N  in 


OI  VO 

d  d 
-H  -H 


in  ™ ~ 

^  in  m 
ra  oj  »n 


oj 

<N 

-H 


o 

-H 


-H 

VO 


-H 


m  vo  in  vo 


d 

-H 


-H 


o 

-H 


-H 

o 


o 

-H 


<u 

s 

-C 

d 

< 

C3 

Oh 

3 

T3 

23 

13 

1/3 

io 

.d 

C/3 

u 

13 

Oh 

JO 

<u 

d 

<u 

d 

u 

o 

CL 

> 

X 

X 

X 

X 

X 

u 

d 

W 

£ 

H 

u 

> 

Q 

X 

X 

X 

<D 

-G 

C* 

X 

X 

X 

K 

S 

X 

X 

X 

X 

> 

^  Springer 


VI  Titan  2574.73  ±  0.09  2575.15  ±  0.02  2574.78  ±  0.06  2574.47  ±  0.06  0.26 


Table  5  continued 


126 


B.  A.  Archinal  et  al. 


&  g 

<D  M 
TO  w 


£ 

M 


I  | 


OC  TO 


23  £ 

o  M 

Oh  w 


'■an  M 

°  -e  c- 

b£)  2  0/5 
3  3 

©  3  ^ 
-,<1  ”  c3 

<  O  2 


>» 

«3  'g' 

U  (1  i 
3  -c  22 

23  o  ^ 

O.  C3  2 

cn  §?  2 


£ 

i2 


in  NO 

-H  -H 

t"; 

CN  CN 


in  os 
TT  CN 

-H  -H 


>  > 


-H  -H 


oo 

CN  © 

-H  -H 

m  in 


3 

<D 

PH 


os  so  os  so  co  *n 


cNcNco^cnin^o^^in 

-H-H-H-H-H-H-H-H-H-H 

min-H-Hr-cN^tr-cnr- 
— 'CNcncN^i-incNcnt^ 


-H  -H 

os  <u  <u 

cn  i>  £  £ 

cn  t--  3  c3 

cn  in  cn  cn 


o 

-H 


-H 

n- 


o  o 
-H  -H 

OO  Os 
in  oo 


o 

U 


3 

2 

O 


CN  ON  1)  1) 


£  £  £ 


cn  in  oo  oo  oo 


£  £ 

cti  c3 

cn 


so  OO  OO  SO 


CN 

-H 

N- 


.2 

o3 

c3 

3 

3 

o 

TD 

3 

3 

TO 

3 

TD 

13 

O 

3 

3 

a 

3 

« 

.3 

K/t 

13 

.2 

3 

0 

-a 

o 

2 

s 

<L) 

Vh 

u 

o 

T3 

<D 

Q 

o 

>— 5 

2 

o 

Oh 

o 

04 

13 

03 

> 

o 

3 

CLh 

> 

C 

3 

— 

s 

13 

"2 

£ 

D 

’£ 

2 

P 

0 

hO 

0 

> 

> 

>< 

X 

X 

>< 

>< 

X 

X 

> 

a 

5 

> 

^  Springer 


Table  5  continued 


Report  of  the  IAU  Working  Group 


127 


&  g 

<D  M 
TD  w 


s 

M 


04  "O 


<2  g 
O  ^ 
Oh  W 


-  s 
.2  ^ 
S-H  ' - ' 

O  W2 
eS  2 


>» 

2  ^  S' 

u  cd  v 

c  -c  w 
2  O  a> 
g-  2  2 

5  ?  2 


VO  OO 

-H  -H  -H  -H 


Oh 

<U 

z 


-H 

VO 


rd 

H 

Oh 

CD 

Q 

rd 

o 

H 

rd 

t-l 

Oh 

> 

> 

> 

> 

> 

^  Springer 


128 


B.  A.  Archinal  et  al. 


With  the  introduction  of  large  mass  storage  to  computer  systems,  digital  cartography 
has  become  increasingly  popular.  Cartographic  databases  are  important  when  considering 
irregularly  shaped  bodies  and  other  bodies,  where  the  surface  can  be  described  by  a  file  con¬ 
taining  the  coordinates  for  each  pixel.  In  this  case  the  reference  sphere  has  shrunk  to  a  unit 
sphere.  Other  parameters  such  as  brightness,  gravity,  etc.,  if  known,  can  be  associated  with 
each  pixel.  With  proper  programming,  pictorial  and  projected  views  of  the  body  can  then  be 
displayed. 

Taking  all  of  this  into  account,  our  recommendation  is  that  longitudes  on  dwarf  plan¬ 
ets,  minor  planets,  their  satellites,  and  comets  should  be  measured  positively  from  0  to  360 
degrees  using  a  right-hand  system  from  a  designated  prime  meridian.  The  origin  is  the  center 
of  mass,  to  the  extent  known. 

Latitude  is  measured  positive  and  negative  from  the  equator;  latitudes  toward  the  positive 
pole  are  designated  as  positive.  For  regular  shaped  bodies  the  cartographic  latitude  of  a 
point  on  the  reference  surface  is  the  angle  between  the  equatorial  plane  and  the  normal  to 
the  reference  surface  at  the  point.  In  the  cartographic  system,  the  position  of  a  point  (P) 
not  on  the  reference  surface  is  specified  by  the  cartographic  latitude  of  the  point  ( P')  on 
the  reference  surface  at  which  the  normal  passes  through  P  and  by  the  height  (h)  of  P 
above  P' . 

For  irregular  bodies  orthographic  digital  projections  often  are  adopted  for  cartographic 
portrayal  as  these  preserve  the  irregular  appearance  of  the  body  without  artificial  distortion. 
These  projections  should  also  follow  the  right-hand  rule. 

Table  6  contains  data  on  the  size  and  shape  of  selected  dwarf  planets,  minor  planets,  their 
satellites,  and  comets.  The  first  column  gives  the  effective  radius  of  the  body  and  an  estimate 
of  the  accuracy  of  this  measurement.  This  effective  radius  is  for  a  sphere  of  equivalent  vol¬ 
ume.  The  next  three  columns  give  estimates  of  the  radii  measured  along  the  three  principal 
axes. 

The  uncertainties  in  the  values  for  the  radii  in  Table  6  are  generally  those  given  by  the 
authors,  and,  as  such,  frequently  have  different  meanings.  Sometimes  they  are  standard  errors 
of  a  particular  data  set,  sometimes  simply  an  estimate  or  expression  of  confidence. 

A  dwarf  planet,  minor  planet,  one  of  their  satellites,  or  a  comet  will  typically  be  included 
in  Table  6  only  if  it  meets  fundamental  publication,  data  quality,  and  applicability  criteria. 
Estimated  values  for  the  body’s  size  and  shape  (modeled  as  a  spheroid)  should  have  been 
published  in  a  refereed  journal  (or  a  revision  that  is  in  preparation).  The  analysis  to  determine 
these  values  will  have  been  derived  from  data  of  sufficient  fidelity  and  quality  to  assure  an 
accurate  estimate,  and  some  portion  of  that  data  will  have  been  acquired  via  direct  methods 
(e.g.  direct  imaging  from  a  spacecraft,  a  space  telescope,  or  an  adaptive  optics  system).  Lastly, 
a  cartographic  need  must  exist  that  justifies  the  definition  of  a  size  and  shape  for  this  body. 

The  radii  given  in  the  tables  are  not  necessarily  the  appropriate  values  to  be  used  in  dynam¬ 
ical  studies;  the  radius  actually  used  to  derive  a  value  for  the  dynamical  form  factor  (J2)  (for 
example)  should  always  be  used  in  conjunction  with  it. 


8  Recommendations 

Although  we  have  not  done  so  in  the  past,  the  Working  Group  briefly  summarizes  here 
what  we  have  come  to  see  as  current  urgent  needs  relative  to  the  development  of  planetary 
cartographic  products. 
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Table  6  Size  and  shape  parameters  of  selected  dwarf  planets,  minor  planets,  their  satellites,  and  comets 


Body 

Effective  radius  ( km) 

Radii  measured  along  principal 
( km)  ( km) 

1  axes 
(km) 

(1)  Ceres 

476.2  ±  1.7 

487.3  ±  1.8 

Same 

454.7  ±  1.6 

(a) 

(4)  Vesta 

289  ±5 

280  ±  5 

229  ±5 

(21)  Lutetia 

52.5  ±  2.5 

62.0  ±  2.5 

50.5  ±  2.0 

46.5  ±  6.5 

(243) Ida 

15.65  ±  0.6 

26.8 

12.0 

7.6 

(253)  Mathilde 

26.5  ±  1.3 

33 

24 

23 

(433)  Eros 

8.45  ±  0.02 

17.0 

5.5 

5.5 

(511)  Davida 

150 

180 

147 

127 

(b) 

(951)  Gaspra 

6.1  ±0.4 

9.1 

5.2 

4.4 

(2867)  Steins 

2.70 

3.24 

2.73 

2.04 

(4179)  Toutatis 

2.13 

1.015 

0.85 

(25143)  Itokawa 

0.535 

0.294 

0.209 

(134340)  Pluto 

1195  +  5 

Same 

Same 

Same 

(134340)  Pluto:  I  Charon 

605  ±  8 

1  P/Halley 

8.0  ±  0.5 

4.0  ±  0.25 

4.0  ±  0.25 

9P/Tempel  1 

3.0  ±0.1 

3.7 

2.5 

(c) 

19P/Borrelly 

4.22  ±  0.05 

3.5  ±  0.2 

- 

- 

81  P/Wild  2 

1.975 

2.7 

1.9 

1.5 

(a)  An  oblate  spheroid.  Carry  et  al.  (2008)  cite  values  of  r  =  467.6  ±  2.2,  a  =  b  =  479.7  ±  2.3,  c  =  444.4  ± 
2.1.  However  the  amount  of  improvement  over  the  Thomas  et  al.  (2005)  values  is  not  clear,  so  no  change  is 
recommended 


(b)  Values  have  been  updated  from  the  diameter  values  which  appear  in  Conrad  et  al.  (2007)  to  the  radii  values 
given  above,  which  result  from  additional  observations  and  appear  in  a  publication  by  the  same  authors  (in 
preparation) 

(c)  The  maximum  and  minimum  radii  are  not  properly  the  values  of  the  principal  semi-axes,  they  are  half  the 
maximum  and  minimum  values  of  the  diameter.  Due  to  the  large  deviations  from  a  simple  ellipsoid,  they  may 
not  correspond  with  measurements  along  the  principal  axes,  or  be  orthogonal  to  each  other 


1.  The  importance  of  geodetically  controlled  cartographic  products — i.e.,  derived  from 
least  squares  photogrammetric,  radargrammetric,  or  altimetric  (cross-over)  solutions — 
is  well  known.  These  products  are  valuable  since  they  are  precise  and  cosmetically  ideal 
products  at  the  sub-pixel  level  of  the  data,  with  known  or  derivable  levels  of  precision 
and  accuracy.  In  addition  global  control  solutions  also  provide  for  improved  body  pole 
position,  spin,  and  shape  information,  with  reduced  effects  of  random  error  and  often 
systematic  error.  Such  solutions  would  allow  for  improvements  in  the  recommended 
models,  and  more  importantly  provide  for  higher  (and  known)  precision  and  accuracy 
cartographic  products.  Although  a  flood  of  new  planetary  datasets  is  currently  arriving, 
it  appears  that  the  production  of  such  products  is  often  not  planned  for  or  funded.  We 
strongly  recommend  that  this  trend  be  reversed  and  that  such  products  be  planned  for 
and  made  as  part  of  the  normal  mission  operations  and  data  analysis  process. 

2.  As  indicated  above,  Konopliv  et  al.  (2006)  have  developed  an  improved  model  for  the 
orientation  of  Mars.  We  note  that  by  the  time  of  our  next  report,  the  orientation  of  Mars 
from  that  model  will  start  to  diverge  significantly  from  that  currently  recommended  by 
this  Working  Group  and  dating  to  our  2000  report.  We  recommend  their  model  be  updated 
or  a  similar  model  be  developed  that  takes  advantage  of  the  substantial  additional  Mars 
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data  available  since  the  time  of  their  work,  so  that  it  can  be  adopted  by  this  Working 
Group  and  operational  Mars  missions  in  the  2012  time  frame. 

3.  As  described  above,  there  appear  to  be  a  number  of  slightly  conflicting  determinations 
for  the  rotation  rates  of  Jupiter  and  Saturn.  We  urge  the  planetary  community  to  jointly 
address  resolving  the  various  determinations  and  to  develop  consensus  determinations, 
such  as  was  done  in  the  past  for  Jupiter  by  Riddle  and  Warwick  (1976). 
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Appendix:  changes  since  the  last  report 

This  appendix  summarizes  the  changes  that  have  been  made  since  the  2006  report  ( Celestial 
Mechanics  and  Dynamical  Astronomy  98,  155-180,  2007). 

1.  Following  the  list  of  dwarf  planets  defined  by  the  IAU  WGPSN  and  the  CSBN  (2010), 
the  entries  for  Pluto  and  Charon  have  been  moved  to  Tables  3  and  6  and  adapted  to  the 
positive  pole  definition,  and  the  entry  for  Ceres  moved  to  Table  6. 

2.  A  new  pole  and  rotation  model  has  been  implemented  for  Mercury,  based  on  the  deter¬ 
mination  of  Margot  2009,  but  using  the  previous  value  of  Wo-  See  Sect.  4  for  further 
details. 

3.  An  algorithm  is  described  in  the  text  for  expressing  the  orientation  of  the  Moon  using 
the  JPL  DE  421  lunar  ephemeris,  rotated  to  the  mean  Earth/polar  axis  system,  in  order 
to  obtain  the  pole  and  rotation  with  high  precision. 

4.  The  new  pole  and  rotation  model  of  Konopliv  et  al.  (2006)  for  Mars  continues  to  be 
noted  in  the  text,  but  is  not  recommended  for  general  use  at  this  time. 

5.  The  rotation  rate  of  Jupiter  was  updated  in  the  2000  report  (Seidelmann  et  al.  2002) 
based  on  a  new  determination  by  Higgins  et  al.  (1997).  That  rate  is  different  from 
the  previously  accepted  value  (“System  III  (1965)”)  of  Riddle  and  Warwick  (1976) 
and  is  not  in  agreement  with  new  determinations  by  Russell  et  al.  (2001)  and  Yu  and 
Russell  (2009).  After  discussion  with  some  of  these  authors  (Higgins  2010,  private 
communication;  Russell  2010,  private  communication),  here  it  is  changed  back  to  the 
previous  Riddle  and  Warwick  (1976)  value.  Higgins,  Russell,  and  the  Working  Group 
recommend  that  a  community  consensus  on  a  new  value  be  achieved — much  as  was 
facilitated  by  Riddle  and  Warwick — before  a  new  value  is  adopted. 

6.  Several  proposed  new  rotation  models  for  Saturn  are  noted  in  the  text,  but  not  recom¬ 
mended  for  general  use  at  this  time.  As  with  Jupiter,  the  Working  Group  recommends 
that  a  community  consensus  on  a  new  model  be  achieved  before  a  new  value  is  adopted. 

7.  The  rotation  models  of  Mimas,  Enceladus,  Tethys,  Dione,  and  Iapetus  have  been 
updated  in  Table  2  based  on  the  Table  24.4  values  in  Roatsch  et  al.  (2009a).  (The 
Dione  model  is  originally  from  Roatsch  et  al.  (2008b);  the  Mimas,  Tethys,  and  Iapetus 
models  are  originally  from  Roatsch  et  al.  (2009b);  the  revised  Enceladus  model  was  not 
yet  derived  when  Roatsch  et  al.  (2008a)  was  published.).  These  changes  are  to  the  values 
of  Wo  for  these  bodies,  based  on  semi-controlled  (including  limb  fitting)  new  mapping 
of  these  bodies.  Note  that  in  some  further  interim  revisions  of  the  relevant  mapping 
products.  Wo  has  not  been  further  updated,  resulting  in  some  variation  in  longitude  of 
the  defining  features  in  the  interim  products.  Final  products  will  be  made  later  with 
associated  further  updates  to  the  Wo  values  (Roatsch  2009,  private  communication). 
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The  rotation  model  of  Titan  has  been  updated  in  Table  2  based  on  the  improved  pole 
position  and  synchronous  spin  model  of  Stiles  et  al.  (2008)  and  as  formulated  in  Cassini 
Project  (2009).  This  model  was  derived  from  Cassini  RADAR  observations  of  Titan. 
A  model  for  nonsynchronous  spin,  resulting  from  surface-atmosphere  interactions  and 
possibly  indicating,  by  its  magnitude,  the  presence  of  a  sub-surface  ocean  (Lorenz  et  al. 
2008),  is  as  yet  not  sufficiently  definitive  to  be  adopted.  Note  that  the  adopted  model 
does  not  include  any  terms  for  precession  and  nutation,  but  for  pole  position  fits  the 
data  better  than  the  previous  WG  model  that  assumed  a  pole  perpendicular  to  the  orbital 
plane.  Stiles  et  al.  (2010)  note  a  correction  to  their  earlier  work,  substantially  reducing 
the  inferred  nonsynchronous  rotation  rate  and  giving  a  slightly  different  pole  position, 
but  they  do  not  provide  a  compatible  spin  model,  so  the  newer  determination  is  not 
adopted  here. 

8.  Clarification  has  been  added  to  Sect.  VI  to  explain  when  positive  longitudes  are  east  or 
west  in  the  planetocentric  and  planetographic  coordinate  systems. 

9.  (1)  Ceres  has  been  added  to  Table  3  based  on  work  described  by  Thomas  et  al.  (2005), 
Li  et  al.  (2006a),  and  Li  et  al.  (2006b).  Wo  has  been  derived  from  W  =  333.14°  at  the 
UTC  epoch  2003-12-28T10:5 1 :59  (Thomas  2010,  private  communication).  The  period 
of  9.074170  hours  from  Chamberlain  et  al.  (2007)  has  been  used  to  derive  the  rotation 
rate.  (2)  Pallas,  (21)  Lutetia,  (5 1 1)  Davida,  and  (2867)  Steins  have  been  added  to  Table  3 
from  Carry  et  al.  (2010a)  (Pallas);  Carry  et  al.  (2010b)  (Lutetia),  Conrad  et  al.  (2007) 
(Davida);  and  Jorda  et  al.  (2010)  and  Keller  et  al.  (2010)  (Steins).  Schmidt  et  al.  (2009) 
have  determined  a  pole  position  for  (2)  Pallas  of  <xq  =  42°  ±  10  and  Sq  =  — 12°  ±  10. 
However,  since  they  provide  no  corresponding  information  on  spin,  values  for  the  pole 
position  and  spin  from  Carry  et  al.  (2010a)  have  been  recommended  here.  The  pole 
position  for  (4)  Vesta  has  been  updated  from  Li  et  al.  (2010). 

10.  The  equatorial  radius  of  the  Sun  has  been  added  to  Table  4  for  comparison,  as  speci¬ 
fied  in  the  IAU  (1976)  System  of  Astronomical  Constants  (International  Astronomical 
Union  (IAU)  1977,  pp.  31,  52-66). 

11.  The  size  and  shape  of  the  Earth  has  been  updated  based  on  International  Association 
of  Geodesy  recommendations  (Groten  2000).  The  (zero  frequency  tide)  mean  equato¬ 
rial  radius,  a,  is  used  as  is  (ibid.,  equation  46),  while  the  (zero  frequency  tide)  mean 
polar  flattening,  1//  (ibid.,  equation  22),  is  used  to  calculate  the  mean  polar  radius, 
b  =  a  —  af .  The  mean  radius  is  calculated  via  the  arithmetic  mean  (2a  +  b) /3  (Moritz 
1980,  392).  The  RMS  deviation  from  the  spheroid,  the  maximum  elevation,  and  max¬ 
imum  depression  has  been  left  at  their  previous  values  since  they  are  only  intended  to 
be  approximate.  The  a  and  1  //  values  used  are  also  those  recommended  in  the  IERS 
Conventions  (McCarthy  and  Petit  2004). 

12.  Sizes  and  shapes  have  been  updated  in  Table  5  for  Phobos  and  Deimos  based  on  Willner 
et  al.  (2010)  and  Thomas  (1993)  respectively. 

13.  Sizes  and  shapes  have  been  updated  in  Table  5  for  the  Galilean  satellites  of  Jupiter  based 
on  Thomas  et  al.  (1998)  (Io);  Anderson  et  al.  (2001a)  (Ganymede),  Anderson  et  al. 
(2001b)  (Callisto),  and  Nimmo  et  al.  (2007)  (Europa).  While  making  these  updates,  we 
have  noted  that  the  previous  recommended  sizes  and  shapes  for  the  Galilean  satellites 
do  not  originate  from  Davies  et  al.  (1998)  as  indicated  in  the  2000  report  (Seidelmann 
et  al.  2002)  but  rather  from  the  Galileo  SSI  instrument  team  (K.  P.  Klassen,  e-mail  of 
2000  November  13).  The  mean  radii  for  Europa,  Ganymede,  and  Callisto  there  came 
from  M.  Davies  and  T.  Colvin  via  photogrammetric  solutions  and  the  mean  radius  for 
Io  and  the  triaxial  radii  for  all  four  satellites  came  from  P.  Thomas  via  limb  fit  solutions. 
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14.  Sizes  and  shapes  have  been  updated  in  or  added  to  Table  5  for  the  Saturnian  satellites 
(in  order  outward  from  Saturn)  Pan,  Daphnis,  Atlas,  Prometheus,  Pandora,  Epimetheus, 
Janus,  Mimas,  Methone,  Pallene,  Tethys,  Telesto,  Calypso,  Dione,  Helene,  Polydeuces, 
Rhea,  Hyperion,  Iapetus,  and  Phoebe,  based  on  new  determinations  listed  by  Thomas 
(2010);  and  for  Anthe  based  on  a  new  determination  listed  by  Porco  et  al.  (2007).  Size 
and  shape  information  for  Titan  is  updated  from  Zebker  et  al.  (2009).  Also  note  that  in 
the  2006  report  (Seidelmann  et  al.  2007,  pp.  174-175)  the  along  orbit  equatorial  radius 
of  Metis  should  have  been  listed  as  20  km,  the  polar  radius  as  17  km,  and  the  RMS 
deviation  from  the  ellipsoid  as  blank;  and  the  mean  radius  of  Helene  should  have  been 
listed  as  16km  and  the  RMS  deviation  from  the  ellipsoid  as  0.7 km. 

15.  In  Table  6  a  mean  radius  value  from  Thomas  et  al.  (2005)  has  been  added  for  ( 1 )  Ceres 
and  a  footnote  added  to  note  an  alternate  determination  of  Ceres’  shape  by  Carry  et  al. 
(2008).  Clarifying  what  was  shown  in  the  2006  report  (Seidelmann  et  al.  2007,  p.  177) 
since  it  is  an  oblate  spheroid,  the  first  and  second  axes  are  listed  as  being  the  same, 
different  from  the  third  axis.  (21)  Lutetia,  (511)  Davida  and  (2867)  Steins  have  been 
added  to  Table  6  from  Drummond  et  al.  (2010)  (Lutetia),  Conrad  et  al.  (2007)  (Dav¬ 
ida),  and  Jorda  et  al.  (2010)  and  Keller  et  al.  (2010)  (Steins).  Schmidt  et  al.  (2009)  and 
Carry  et  al.  (2010a)  have  both  determined  shape  parameters  for  (2)  Pallas.  These  are 
respectively:  r  =  272,  a  =  291  ±  9,  b  =  278  ±  9,  and  c  =  250  ±  9;  r  =  256  ±  3,  a  =  275 
±  4,  b  =  258  ±  3,  and  c  =  238  ±  3  km.  We  have  been  unable  to  clearly  differentiate 
which  is  superior  so  make  no  recommendation  as  to  which  to  use  at  this  time. 
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